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rho  bistatic  radar  cross  section  of  a  perfectly  conducting  flat  plate  that  has 
a  rhombus  shape  (equilateral  parallelogram)  is  investigated.  The  Ohio  State  Uni¬ 
versity  electromagnetic  surface  patch  code  (ESP  version  4)  is  used  to  compute  the 
theoretical  bistatic  radar  cross  section  of  a  35-  x  27-in  rhombus  plate  at  1.3  GHz 
over  the  bistatic  angles  15°  to  142?.  The  ESP-4  computer  code  is  a  method  of  mo¬ 
ments  FORTRAN-77  program  which  can  analyze  general  configurations  of  plates 
and  wires.  This  code  has  been  installed  and  modified  at  Lincoln  Laboratory  on  a 
SUN  3  computer  network.  Details  of  the  code  modifications  are  described.  Com¬ 
parisons  of  the  method  of  moments  simulations  and  measurements  of  the  rhombus 
plate  are  made.  It  is  shown  that  the  ESP-4  computer  code  provides  a  high  degree 
of  accuracy  in  the  calculation  of  copolarized  and  cross-polarized  bistatic  radar  cross 
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1.  INTRODUCTION 


In  radar  cross  section  (RCS)  measurements  of  complex  targets  it  is  desirable  to  make  additional 
measurements  of  reference  targets  such  as  flat  plates  and  cylinders  for  which  the  radar  cross  section 
is  readily  computed.  These  reference  targets  are  used  in  checking  the  calibration,  target  quiet-zone 
characteristics,  and  mechanical  alignment  of  the  measurements  system.  This  report  addresses  the 
simulated  RCS  of  a  perfectly  conducting  flat  plate  used  in  bistatic  measurements.  Depicted  in 
Figure  1-1  is  a  typical  geometry  for  bistatic  RCS  measurements.  The  angle  0  denotes  the  bistatic 
angle  which  is  fixed,  and  the  angle  0  denotes  the  target  rotation  angle  or  aspect  angle. 

•4?Br  ■ 


Figure  1-1.  Geometry  for  bistatic  radar  cross  section  measurements.  The  angle  0  is 
the  bistatic  angle,  and  the  angle  9  is  the  target  rotation  angle. 

Except  for  spheres  and  cylinders,  little  published  data  exist  for  the  far-field  bistatic  radar 
cross  section  of  targets  [1-3].  The  bistatic  scattering  from  a  sphere  is  readily  computed  [4]  and  is 
often  used  as  a  primary  calibration  standard  in  RCS  measurements.  Bistatic  scattering  analyses 
and  measurements  for  a  conducting  cylinder  of  finite  length  are  described  in  [5*8].  Monostatic 
scattering  from  square  and  circular  flate  plates  is  treated  thoroughly  in  [1-3].  A  physical  optics 
formulation  for  the  bistatic  scattering  of  a  polygonal  flat  plate  is  reviewed  in  [3] .  The  bistatic  RCS 
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of  rectangular  and  triangular  plates  has  been  calculated  in  [8].  The  purpose  of  this  report  is  to 
compare  simulated  and  measured  data  for  a  rhombus-shaped  (equilateral  parallelogram)  flat  plate 
over  a  broad  range  of  bistatic  angles. 

A  method  of  moments  code  called  ESP-4  (electromagnetic  surface  patch  code:  version  4)  which 
is  capable  of  analyzing  a  wide  variety  of  antenna  or  radar  cross  section  problems  has  been  devel¬ 
oped  at  The  Ohio  State  University  [9).  The  software  can  analyze  complex  geometries  involving 
multipie-connected  plates  and/or  wires  and  is  well-suited  to  analyzing  the  scattering  from  an  iso¬ 
lated  quadrilateral  plate  like  the  rhombus.  Many  of  the  subroutines  in  the  ESP-4  code  are  based  on 
subroutines  from  an  earlier  wire-grid  RCS/antenna  method  of  moments  code  developed  by  Rich¬ 
mond  [10].  Wire-grid  modeling  of  a  continuous  surface,  such  as  a  fiat  plate,  is  approximated  with 
closely  spaced  conducting  wires  [11].  A  similar  wire-grid  moment  method  computer  code  is  de¬ 
scribed  in  [12].  In  contrast,  surface  patch  modeling  allows  a  piecewise-continuous  approximation  to 
a  complicated  surface.  Surface  patch  modeling  requires  fewer  unknowns  than  wire-grid  modeling: 
hence,  larger  surfaces  can  be  modeled.  ESP-4  utilizes  the  electric  field  integral  equation  (EFIE) 
to  enforce  the  boundary  condition  of  the  tangential  electric  field  being  zero  at  the  surface  of  the 
antenna/target  of  interest.  The  EFIE  solution  allows  lor  either  open  or  closed  surfaces.  The  basis 
and  testing  functions  used  in  this  code  are  piecewise  sinusoidal.  The  surface  patches  (or  surface 
current  modes)  are  assumed  to  have  zero  thickness  and  are  quadrilateral  dipoles  which  are  useful 
in  modeling  an  arbitrary-shaped  body.  For  a  complete  description  of  the  theory  and  capabilities 
of  the  ESP-4  code,  the  reader  is  referred  to  the  user's  manual  [9].  Other  moment  method  codes 
exist  which  are  based  on  the  magnetic  field  integral  equation  (MFIE)  for  closed  surfaces  [13,14] 
and  EFIE  for  arbitrary  surfaces  [15-17].  The  EFIE  moment  method  formulation  used  in  [15-17] 
differs  from  the  ESP-4  formulation  in  that  triangular-shaped  basis  functions  with  pulse  weighting 
are  used  for  the  surface  current  mode  vectors. 

This  report  is  organized  in  the  following  manner:  In  Section  2,  details  of  the  modifications 
made  to  the  electromagnetic  surface  patch  code  ESP-4  are  given.  The  basic  modifications  of  the 
ESP-4  code  are  a  change  in  the  input  data  structure  and  a  new  option  for  fixed  bistatic  angle  R.CS 
patterns.  A  listing  of  the  revised  ES'’-4  main  program  is  given  in  the  appendix.  Bistatic  RCS 
patterns  for  a  35-  x  27-in  rhombus  flat  plate  at  1.3  GHz  are  computed  with  the  ESP-4  code  and 
the  results  are  given  in  Section  3.  Bistatic  angles  from  15°  to  142°  have  been  considered  as  well 
as  two  plate  orientations  (untilted  and  tilted).  The  method  of  moments  simulations  are  compared 
against  measured  data  and  good  agreement  is  demonstrated. 
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2.  ELECTROMAGNETIC  SURFACE  PATCH  CODE  (ESP-4) 

MODIFICATIONS 


The  purpose  of  this  section  is  to  describe  the  modifications  made  to  the  main  program  of 
The  Ohio  State  University  (OSU)  electromagnetic  surface  patch  code  (ESP-4).  This  software  was 
obtained  in  July  1988  from  OSU  and  has  since  been  modified  at  Lincoln  Laboratory.  The  important 
modifications  are  the  addition  of  namelist  input  data  and  an  option  for  bistatic  radar  cross  section 
calculation  with  a  fixed  bistatic  angle  and  variable  target  rotation.  A  complete  description  of  how 
to  use  the  ESP-4  code  can  be  found  in  the  user's  manual  [9].  The  revisions  to  the  ESP-4  code  are 
summarized  graphically  in  Figure  2-1  and  are  described  in  detail  in  Sections  2.1  and  2.2. 
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Figure  2-1.  Block  diagram  depicting  the  revisions  made  to  the  ESP-4  computer  pro¬ 
gram.  The  input  data  have  been  restructured,  and  a  new  option  for  fixed  bistatic  angle 
RCS  pattern  computation  has  been  added. 


2.1  MODIFICATIONS  TO  INPUT  DATA 


The  original  ESP-4  code  defined  plate/'wire  geometries  and  antenna/RCS  pattern  cuts  through 
a  series  of  free- format  |REAI)(  1 1  .*)]  read  statements  on  device  11.  There  are  nearly  80  input 
parameters  in  the  input  data  file  and  the  addition  of  FORTRAN  NAMELISTS,  where  every  input 
quantity  is  clearly  defined,  is  very  attractive.  When  a  NAMELIST  is  used,  the  input  data  are  in 
free  format  and  a  typical  input  parameter,  say  A,  is  of  the  explicit  form  A=1.0,  for  example.  The 
NAMELIST  data  file  is  helpful  in  constructing  input  files  for  the  ESP-4  program  that  are  user- 
friendly  and  easily  debugged.  The  present  Lincoln  Laboratory  version  of  ESP-4  has  attempted  to 
keep  the  definition  of  the  input  paiameters  the  same  as  in  the  OSU  ESP-4  user's  manual  [9] .  The 
manual  should  be  consulted  for  the  definition  of  most  of  the  input  parameters.  Here,  only  newly 
defined  parameters  (input  variables)  will  be  defined.  Nine  namelists  have  been  added  to  the  ESP-4 
code  and  all  READ(1 1.*)  statements  have  been  commented,  The  NAMELIST  reads  have  the  form 
READ(ll.NAM)  where  NAM  is  the  name  of  the  NAMELIST.  The  added  namelist  statements  (in 
FORTRAN)  are  summarized  below  and  follow  line  ESP01G50  of  the  ESP-4  code. 

NAMELIST  /RNCTRL/NGO , NPRINT , NRUNS , NWGS ,  I VR , IWRZT , I NT , INTP , INTO , 
2INWR,IRGM,IFIL,RF,INDZI 

NAMELIST  /FSWEEP/FMC1 ,FMC2 , DFZI ,DFF , IRS 12 ,THRD ,PHRD ,THRI ,PHRI 
NAMELIST  /PATTRN/IFE , 1PFE.KNDFE , PHFE.IFA . IPFA , FNDFA , THFA , 

2ISE , IPSE .FNDSE , PHSE ; THIN , PHIN , ISA , IPSA .FNDSA , THSA , 

3AZRANG , ELRANG , AZMIN , ELM IN .IBISC, BETA , NPTBIS , BANGRG 
NAMELIST  /FWIRET/FMC , CMM , A , NPLTS 

NAMELIST  /PLATEG/NCNRS , SEGK , IREC , IPN , IGS , ZSHT , XP , YP , ZP 
NAMELIST  /SAVEZ/IWRZM, IRDZM 

NAMELIST  /WIREAG/NM , NP , N AT , NFPT , NFS1 ,NFS2 , X , Y , Z , I A ,  IB 

NAMELIST  /GENLOD/IFMM . I ABB . VLGG , ZLL 

NAMELIST  i ATTACH/NASAT , I ABAT , NPLA , VGA , ZLDA ,BDSK 

The  relation  between  the  new  namelists  and  the  old  read  statements  is  as  follows: 

•  RNCTRL=READ  #1 

«  FSWEF.P  =  READ  #1A 

•  PATTRN=READS  #  2.3,4.5  with  the  addition  of  parameters  AZRANG,  EL¬ 
RANG.  AZMIN.  ELMIN.  IBISC,  BETA.  NPTBIS.  BANGRG  described  in  Sec¬ 
tion  2.2. 

•  FWIRET=READ  #6  with  the  addition  of  the  parameter  NPLTS  (number  of 
plates)  from  READ  #7.  Namelist  PLATEG  (described  next)  is  read  NPLTS 
number  of  times. 

•  PLATEG=READS  #  7,  8,  9  with  the  substitution  of  explicit  (XP,  YP,  ZP) 
coordinates  in  array  format  for  the  corners  of  plate  NPL.  The  software  then 
fills  in  the  original  three-dimensional  PCN  arrays  (see  ESP-4  manual)  which 
contain  the  plate  corner  coordinates. 
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•  SA\VEZ=READ  #10 

•  WIREAG- READS  #11.  12.  13 

•  GENLOD=READ  #14  is  executed  NFPT  times  corresponding  to  the  number 
of  feed  points. 

•  ATTACH=READ  #15  is  executed  NAT  times  corresponding  to  the  number  of 
attachment  points  between  plates  and  wires, 


2.2  MODIFICATIONS  FOR  FIXED  BISTATIC  ANGLE  RCS  PATTERNS 

Consider  Figure  2-2  which  shows  the  bistatic  geometry  for  an  arbitrary  target.  The  pair  of 
angles  (0,.<z>,-)  denotes  the  incident  direction  and.  similarly,  [9a.&a)  denote  the  scattering  direction. 
The  original  version  of  ESP-4  assumed  in  a  bistatic  calculation  that  the  angle  of  incidence  was  fixed 


Figure  2-2.  Geometry  for  bistatic  radar  cross  section  computation  in  the  ESP-4  com¬ 
puter  program. 


and  the  scattering  direction  (or  observation  angle)  was  variable.  Thus,  the  software  computes  the 
bistatic  RCS  pattern  as  a  function  of  bistatic  angle.  As  described  earlier  in  Figure  1-1,  this  differs 
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from  a  typical  bistatic  RCS  measurement  where  the  bistatic  angle  is  fixed  and  the  target  rotation 
angle  varies.  The  Lincoln  Laboratory  version  of  ESP-4  uses  a  new  option,  IBISC=1,  which  fixes 
the  bistatic  angle  BETA  and  increments  the  incident  and  observation  angles  appropriately.  The 
implementation  of  the  new  option  occurs,  primarily,  following  line  ESP14200.  Note:  If  IB1SC-0  or 
is  not  specified,  then  the  code  performs  in  the  original  bistatic  mode  where  the  angle  of  incidence  is 
fixed  and  the  observation  angle  varies.  When  IBISC=1,  NPTBIS  is  the  number  of  target  rotation 
angles  and  BANGRG  is  the  angular  range  of  target  rotation  in  degrees,  BANGRG  can  be  a 
positive  or  negative  number,  as  appropriate,  to  determine  the  target  rotation  direction  (clockwise  or 
counterclockwise),  The  original  scattering  indicators  ISE,  ISA  must  be  set  equal  to  the  appropriate 
value  to  invoke  bistatic  RCS  computation  in  the  appropriate  elevation  or  azimuth  plane.  For 
example,  for  scattering  in  the  elevation  plane  (<t>  is  constant)  the  following  three  linos  of  input  data 
are  appropriate: 

ISE-2 , IPSEal , FNDSE=3 , 0 ,PHSEa0 . 0 .THIN-90 . 0 ,PHINa0 . 0 , ELfUNGa360 . , 

ISA=0 , IPSAsl , FNDSAs3 . 0 ,THSAa90 . 0 , 

IBISC*=1  .BETA-120.  ,NPTBIS  =  121,BANGRGa~360.  , 

ISE  is  set  equal  to  2  to  indicate  that  bistatic  scattering  in  the  elevation  plane  is  desired.  ISA  is  set 
equal  to  0  to  indicate  that  an  azimuth  plane  scattering  pattern  is  not  desired.  IBISC  is  set  equal 
to  1  to  invoke  the  fixed  bistatic  angle  option  with  BETA==120,  the  bistatic  angle  >n  degrees;  NPT- 
BIS= 12 1 ,  the  number  of  target  rotation  angles  for  the  fixed  bistatic  angle;  and  BANGRG— -360, 
the  angular  range  of  target  rotation  in  the  counterclockwise  direction  (+360  would  produce  clock¬ 
wise  rotation).  With  IEISC=1.  PHIN— 0.0  in  the  first  data  line  is  used;  however,  the  value  of  THIN 
is  computed  within  the  new  section  of  code.  The  initial  value  of  THIN  is  equal  to  -BETA/2  and  is 
then  uniformly  incremented  within  a  new  DO  LOOP  (DO  1920  IBIS-1, NPTBIS)  over  the  target 
rotation  angular  range  BANGRG.  The  initial  observation  angle  or  scattered  angle  is  initially  set 
equal  to  -t-BETA/2  and  is  similarly  incremented.  If  the  target  is  a  flat  plate  located  in  the  xy 
plane,  then  the  first  point  computed  by  the  ESP-4  code  is  the  bistatic  specular  response.  This  fact 
is  due  to  the  broadside  direction  being  chosen  as  the  bisector  of  the  bistatic  angle.  This  choice  is 
arbitrary  and  can  be  changed  as  desired  within  the  revised  ESP-4  code  between  lines  ESP14200 
and  ESP14240.  Note:  Tn  the  above  three  data  lines  if  IBISC  is  set  to  0,  then  the  program  would 
compute  a  bistatic  pattern  in  the  elevation  plane  based  on  the  values  given  in  the  first  data  line 
above.  Thus,  the  program  would  run  as  in  the  original  version  of  ESP-4. 

The  original  ESP-4  code  assumes  that  azimuth  and  elevation  antenna/RCS  patterns  cover  a 
full  360°  angular  range.  New  parameters,  AZRANG,  ELRANG  have  been  included  in  the  input 
data  to  cover  an  arbitrary  azimuth  and  elevation  range,  respectively.  Additionally,  to  start  patterns 
at  an  arbitrary  azimuth  or  elevation  angle,  AZMIN  and  ELMIN,  respectively,  have  been  defined. 
The  default  values  of  AZMIN  and  ELMIN  are  zero  so  that  they  do  not  have  to  be  defined  in  the 
input  data  file.  Further,  the  line  involving  IBISC  does  not  have  to  be  defined  in  the  input  data 
since  appropriate  default  values  are  specified  prior  to  the  READ(ll.PATTRN)  statement.  It  should 
be  noted  that  when  the  new  bistatic  option  is  used,  the  parameters  AZRANG  and  ELRANG  are 
automatically  set  equal  to  zero.  This  effectively  forces  the  number  of  scattering  angles  equal  to  one 
for  a  given  incident  angle.  Other  minor  changes  to  the  code  are  documented  with  comments  in  the 
modified  ESP-4  code  listed  in  the  appendix. 
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3.  RESULTS 


3.1  RHOMBUS  PLATE  MODEL 

For  a  general  rhombus  flat  plate,  as  shown  in  Figure  3-1,  let  L  and  H'  represent  the  diagonal 
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Figure  3-L  Geometry  for  a  rhombus  fiat  plate.  The  rhombus  characteristics  (side 
length,  interior  angles,  and  area]  are  readily  computed  from  the  diagonal  length  L  and 
diagonal  width  IV. 


length  and  diagonal  width,  respectively,  and  let  a  be  the  side  length.  The  side  length  of  the  rhombus 
is  computed  according  to 


VffTTv2 

A  —  . — . - - 

2 

and  the  area  A ,  is  given  by 

A  =  \LW. 

2 


(3.1) 


(3.2) 


The  rhombus  interior  angles  at  and  aw  are  computed  using 
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(3.3) 


-I/* 


ql  =  2sin-1(  — ) 


«iv  =  2 sin  — ). 


(3.4) 


The  specular  monostatic  radar  cross  section  of  the  flat  plate  is  given  by  the  well-known  equation 


a  — 


4tt.42 


(3.5) 


A  sketch  of  the  rhombus  flat  plate  target  under  consideration  is  shown  in  Figure  3-2.  The  plate 
has  a  diagonal  length  L  -  35  in  and  a  diagonal  width  H*  =  2?  in.  The  length  of  each  of  the 


Figure  3-2.  Rhombus  plate  with  diagonal  length  L  -  35  in  and  diagonal  width  W  - 
27  in  used  both  in  ESP-4  simulations  and  in  measurements. 


four  sides  is  approximately  22.1  in.  The  acute  angles  of  this  plate  are  an-  -  75.3°  and  the  obtuse 
angles  are  oi  —  104.7°.  In  terms  of  wavelength,  at  the  frequency  of  interest  1.3  GHz,  the  electrical 
dimensions  of  the  plate  are  L  —  3.85A  and  W  =  2.97A  and  the  plate  area  is  5.74  square  wavelengths. 
The  thickness  of  the  experimental  plate  is  0.5  in  or  0.055A  and  the  plate  is  simulated  using  zero 
thickness.  For  the  given  plate  dimensions,  the  specular  monostatic  radar  cross  section  is  computed 
(using  Equations  3.2  and  3.5)  to  be  13.4  dBsm  at  the  desired  frequency  1.3  GHz.  Figure  3-3  shows 
a  three-dimensional  view  of  the  plate  located  in  the  rectangular  coordinate  system  for  two  plate 
orientations.  In  Figure  3-3(a),  the  plate  is  in  the  xy  plane  and  in  Figure  3-3(b),  the  plate  has  been 
rotated  by  45°  with  respect  to  the  y  axis.  In  this  report,  the  RCS  pattern  cuts  are  always  taken  in 
the  xz  plane.  Notice  that  the  long  dimension  of  the  plate  (35  in)  is  oriented  in  the  x  direction.  A 
two-dimensional  view  (looking  along  the  z-axis)  of  the  plate  in  untilted  and  tilted  configurations  is 
shown  in  Figure  3-4.  The  bistatic  geometry  for  the  plate  is  depicted  in  Figure  3-5.  Note  that  the 
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\*m:  6 


(a) 


(b) 


NOTE:  ALL  DIMENSIONS  ARE  IN  INCHES 

Figure  .3-3.  Three-dimensional  view  of  rhombus  plate  in  (a)  untilted  and  (b)  tilted 
orientation  (rotated  45 4  about  the  x  axis).  All  RCS  patterns  in  this  report  are  taken 
in  the  xz  plane. 


angle  0  represents  the  effective  target  rotation  angle  for  the  fixed  bistatic  angle  0.  When  6  -  0°, 
the  perpendicular  direction  to  the  plate  is  at  the  bisector  of  the  bistatic  angle. 

The  ESP-4  code  was  used  to  analyze  the  bistatic  scattering  patterns  of  the  rhombus  plate. 
Bistatic  angles  of  15,  45,  90,  120,  and  142°  were  considered  corresponding  to  available  35-  x  27-in 
rhombus  flat  plate  measured  RCS  data  collected  in  1987  at  the  US  Air  Force  Radar  Target  Scatter 
Facility  (RATSCAT)[18j.  The  important  input  data  for  the  ESP-4  RCS  simulations  are  listed  below 
with  a  description  of  the  parameters. 
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Figure  3-4.  Two-dimensional  view  (looking  along  the  x-axis)  of  the  plate  in  untilted 
and  tilted  configurations. 


ARNCTRL  (Run  control  parameters) 

NG0*1 ,  NPRINT*2 .  NRUNS= 1 ,  NWGS*  1 ,  IWR-0 ,  IVRZT-0 , I  NT-4 ,  INTP-6 ,  INTD=  18 , 

INWR^O , IRGM= 1 , IFIL=0,RF=- 1.0, INDZI-0 , 

It 

In  NAMELIST  RNCTRL,  the  parameter  NGO-1  runs  the  ESP-4  program  from  start  to  finish. 
Printout  of  input  parameters  and  target  geometry  is  implemented  with  NPRINT=2.  The  program 
is  executed  one  time  (NRUNS=1)  and  NWGS  is  not  used.  The  paramet*  -s  IWR  and  IWRZT  are 
both  set  to  zero  which  means  that  the  induced  modal  currents  and  impedance  matrix  elements, 
respectively,  are  not  printed  out.  The  parameter  INT  is  not  used.  The  number  of  Simpson’s  rule 
integration  intervals  used  in  integrating  over  the  surface  patch  monopoles  is  specified  by  INTP— 6. 
The  parameter  INTD  is  not  used  and  INWR=0  means  that  there  are  no  wires  in  the  target  geometry. 
The  parameter  IRGM  is  not  used.  The  parameter  IFIL-0  means  that  full  surface  patch  test  modes 
are  used,  and  RF=-1  performs  a  far-field  computation.  The  parameter  INDZI  is  set  to  zero  which 
means  do  not  perform  a  frequency  sweep  computation. 

APATTRN  (Pattern  specifications) 

IFE-0 , IPFE* 1 , FNDFE-3 . 0 , PHFE-90 . 0 , 

IFA-0 , IPFA- 1 , FNDFA-3 . 0 . THFA-90 . 0 , 

ISE-2 , IPSE- 1 . FNDSE-3 . 0 , PHSE-0 . 0 , THIN-90 . 0 , PHIN-0 . 0 , ELRANG-0 . , 
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Figure  3-5.  Bistatic  geometry  for  fiat  plate.  The  bistp.tic  angle  is  fixed  and  the  target 
aspect  angle  0  varies.  Note  that  6  =  0:  corresponds  to  the  specular  bistatic  RCS. 


ISA«0 , IPSA= 1 , FNDSA=3 . 0 . THSA=90 . 0 , 

IBISC=1 ,BETAa120. ,NPTBIS=121 ,BANGKGs-360. , 

k 

The  third  and  fifth  lines  of  the  PATTRN  NAMELIST  indicate  that  a  bistatic  scattering  eleva¬ 
tion  pattern  (ISE=2)  with  the  new  fixed  bistatic  angle  option  (IBISC=1)  is  used.  The  bistatic 
angle  is  defined  by  BETA=120  (degrees).  The  azimuth  cut  angle  (<p)  is  PHSE^O  (degrees),  which 
produces  a  pattern  in  the  xz  plane.  The  incident  wavefront  azimuth  angle  (<&)  is  PHIN=0  (de¬ 
grees).  The  number  of  target  rotation  angles  is  NPTBIS=121,  and  the  range  of  target  rotation 
is  BANGRG=-360  (counter  clockwise  rotation).  Lines  1,  2,  and  4  of  PATTRN  indicate  that  no 
other  pattern  cuts  are  desired. 

AFVJIRET  (Frequency,  uire  type,  and  number  of  plates) 

FMC-1300 . 0 , CMM-38 . 0 , A*0 . 00 1 , NPLTS- 1 , 

A 

The  frequency  has  been  set  to  1.3  GHz  (FMC=1300,  frequency  in  MHz)  and  CMM  and  A  are  not 
used.  The  number  of  plates  is  NPLTS— 1. 
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ftPLATEG  (Plate  geometry) 

NCNRS{1)*4 ,SEGM(i)«0.2 ,IREC(l)sO,IPN(l)a3,IGS(l)*0,ZSHT(l)"(0.0,0.0) , 

XP(l)a“. 4446,0, , .4445,0. , 

YP(l)-0. ,-.3429,0. , .3429, 

ZP(l)-0. 0,0. 0,0. 0,0.0, 

ft 

There  are  4  corners  on  the  plate  [NCNRS(1)=4],  and  the  maximum  surface  patch  segment  size  is 
0.2 A.  The  parameter  IREC(1)=0  means  that  the  plate  is  polygonal  rather  than  rectangular.  The 
parameter  IPN(1 )— 3  means  that  two  polarizations  are  used  in  the  surface  patch  currents.  The 
parameter  IGS=0  means  that  ESP-4  selects  a  reference  side  of  the  plate.  A  perfectly  conducting 
plate  is  defined  by  the  complex  sheet  impedance  ZSHT(1)=(0. 0,0.0).  The  four  corners  of  the 
plate  are  described  in  rectangular  coordinates  by  the  data  arrays  XP,  YP,  and  ZP  in  meters  and 
correspond  to  the  coordinates  shown  in  Figure  3-3(a). 


ftSAVEZ  (Save  or  reuse  impedance  matrix) 

IWRZM«1 , IRDZM=0 , 

ft 

The  program  saves  the  impedance  matrix  for  additional  runs  using  the  parameter  IWRZM=1. 

After  performing  an  initial  run  and  storing  the  impedance  matrix,  a  second  run  is  made  in 
which  the  impedance  matrix  is  read-in  using  the  following  input  data: 

ftSAVEZ 

IWRZM=0 , IRDZM=3 , 

ft 

For  the  case  of  the  4oc  tilted  rhombus  plate  [Figure  3-3(b)]  the  input  plate  geometry  NAMELIST 
is  given  by 

ftPLATEG 

NCNRS ( 1 ) =4 , SEGM ( 1 ) =0 . 2 , IREC ( 1 ) a0 , IPN ( 1 ) =3 , IGS ( 1 ) -0 , ZSHT (l)a(0. 0,0.0), 

XP(1)*-. 4445,0. , .4445,0. , 

YP(1)*0 . ,-.2425,0. , .2425, 

ZP(l)a0.0,. 2425, 0.0, -.2425, 

ft 

The  first  step  in  running  the  ESP-4  code  is  to  verify  the  target  (flat  plate)  input  geometry.  This 
is  done  by  specifying  NGO=0  in  the  input  data  file.  After  running  the  ESP-4  program,  an  out¬ 
put  file  is  generated  which  contains  the  input  geometry  information.  Geometry  plotting  software 
called  ESP4GM  was  included  with  the  ESP-4  software  and  was  used  to  generate  visual  information 
about  the  plate  simulation.  ESP4GM  is  a  FORTRAN  77  program  which  utilizes  GKS  (Graphic 
Kernel  System)  software.  The  results  are  summarized  in  Figures  3-6  to  3-9.  Figure  3-6  show's 
the  layout  of  the  overlapping  piecewise-sinusoidal  surface  patch  modes  for  both  polarizations.  For 
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Figure  3-6.  Dual-polarized  surface  patch  layout  generated  by  ESP-4  for  the 
35 -  x  27- in  rhombus  plate  at  1.3  GHz. 


each  polarization  there  are  156  modes  arranged  in  a  12  x  13  grid.  In  terms  of  the  plate  surface 
area  (5.74A2)  there  are  approximately  27  modes  (of  one  polarization)  per  square  wavelength.  Each 
dipole  surface  patch  mode  is  a  parallelogram,  and  two  of  these  modes  arranged  side  by  side  form 
a  rhombus,  which  is  similar  to  the  rhombus  plate.  A  detailed  (enlarged)  view  of  two  contiguous 
patches  of  each  polarization  is  depicted  in  Figure  3-7.  The  polarization  vectors  of  the  two  different 
modes  (which  are  not  quite  orthogonal)  make  an  obtuse  angle  with  respect  to  each  other  that  is 
104.7°  as  determined  from  Equation  (3.4).  One  of  the  dipole  surface  patches  is  made  up  of  two 
monopole  surface  patches,  each  of  which  has  a  rhombus  shape.  In  this  figure,  the  side  of  one  of  the 
rhombus-shaped  monopole  surface  patches  has  a  length  of  0.1 87A.  This  is  found  by  dividing  the 
side  of  the  rhombus  plate,  which  is  2.43A,  into  13  equal  segments.  Note  that  the  computer  program 
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Figure  3-  7.  Enlarged  view  of  contiguous  rhombus-shaped  surface  patch  monopoles 
from  Figure  3-6. 


has  chosen  the  size  of  the  patches  subject  to  the  maximum  patch  size  specified  by  the  parameter 
SEGM(l)-0.2  (wavelengths).  A  three- view  plot  of  the  rhombus  plate  geometry  is  shown  in  Figure 
3-8  for  the  untilted  plate  and  in  Figure  3-9  for  the  45°  tilted  plate. 


3.2  COMPARISON  OF  SIMULATED  AND  MEASURED  BISTATIC  RCS 

Using  the  surface  patch  basis  function  layout  given  in  Figure  3-6,  the  bistatic  radar  cross 
section  of  the  rhombus  plate  at  1.3  GHz  was  computed  by  setting  NGO=l.  An  initial  run  was 
made  where  the  method  of  moments  mutual  impedance  matrix  is  computed  and  stored  in  a  disk 
file  for  later  runs.  A  collection  of  simulated  (with  ESP-4)  and  measured  RCS  patterns  of  the 
rhombus  fiat  plate  covering  bistatic  angles  15°  to  142°  are  shown  in  Figures  3-10  to  3-27.  For  all 
bistatic  angles  both  plate  tilt  angles  (0°  and  45°)  were  used  except  for  the  /?  -  15°  case  where 
measurements  were  available  only  for  the  untilted  case.  For  the  measured  data,  the  notations  H 
and  V  refer  to  horizontal  and  vertical  polarization,  respectively.  A  pair  of  letters,  for  example 
HV,  refers  to  horizontal  transmit  polarization  and  vertical  receive  polarization.  In  the  ESP-4  code, 
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Figure  3-8.  Three-new  plot  of  the  rhombus  plate  in  the  untilted  orientation. 


spheiical  components  are  computed;  the  6  component  is  equivalent  to  H  and  the  <j>  component 
is  equivalent  to  V.  Thus,  the  relation  between  the  four  possible  components  in  the  measured  and 
ESP-4  coordinate  systems  is  (HH,  VV,  HV.  VH)<->  (06,  d>&,  0<t>,  09).  For  the  remainder  of  the  report, 
the  RCS  components  will  be  described  using  the  H  and  V  notations. 

Consider  first  Figures  3-10  and  3-11  which  show  the  RCS  patterns  for  the  four  polarizations 
(HH,  VV,  HV,  VH)  at  the  bistatic  angle  15°.  Notice  in  Figure  3-10  that  the  simulated  copolarized 
(HH,  VV)  components  are  in  good  agreement  with  the  measurements.  The  simulated  specular 
bistatic  RCS  (13.1  dBsm)  agrees  closely  with  the  theoretical  specular  monostatic  RCS  (13.5  dBsm) 
computed  earlier.  In  Figure  3-11,  the  simulated  cross-polarized  RCS  (HV,  VH)  is  below  the  -45- 
dBsm  level  and  so  the  measured  data  are  corrupted  by  the  polarization  isolation  and  background 
clutter  in  the  measurements.  For  all  the  bistatic  angles,  when  the  plate  is  not  tilted  the  simulated 
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Figure  3-9.  Three-view  plot  of  the  rhombus  plate  in  the  45‘  tilted  orientation. 

cross-polarized  RCS  is  very  low,  and  so  the  measurements  are  dominated  by  isolation/clutter- 
induced  errors.  These  cross-polarized  data  appear  in  Figures  3-14,  3-18,  3-22,  and  3-26  and  will 
not  be  discussed  further. 

A  series  of  patterns  for  bistatic  angles  45°  to  142°  are  shown  in  Figures  3-12  to  3-27,  where  both 
plate  tilt  angles  are  considered.  For  example,  Figure  3-12  is  for  the  HH  and  VV  RCS  components 
at  the  45°  bistatic  angle  with  no  target  tilt,  while  Figure  3-13  is  for  a  plate  tilt  angle  of  45°.  Notice 
that  the  effect  of  the  plate  tilt  is  to  break  up  the  main  lobe  in  both  HH  and  VV  components.  Good 
agreement  between  the  measurements  and  simulations  is  evident  for  the  copolarized  components 
for  this  case.  Next,  in  Figures  3-14  and  3-15  the  cross-polarized  RCS  is  shown  for  the  untilted  and 
tilted  plate,  respectively.  The  effect  of  tilting  the  plate  is  to  raise  the  theoretical  cross-polarized 
response;  thus,  the  measured  data  bear  a  resemblance  to  the  simulations  as  in  Figure  3-15.  The 
basic  lobing  structure  is  the  same  for  the  measurements  and  the  simulations,  but  the  amplitude 
of  individual  lobes  differs  by  as  much  as  10  dB.  Next,  for  the  90°  bistatic  angle  case  shown  in 
Figure  3-16  there  is  a  clear  broadening  of  the  main  beamwidth  in  HH  and  VV  polarizations,  and 
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the  simulations  and  measurements  are  in  good  agreement.  Good  agreement  is  also  achieved  in 
Figures  3-17  and  3-19  for  the  tilted  plate,  Figures  3-20  to  3-23  show  the  0  =  120°  case  results. 
The  main  lobe  continues  to  broaden  for  HH  and  VV  components  as  shown  in  Figure  3-20.  There 
is  a  somewhat  better  agreement  between  simulations  and  measurements  for  HH  compared  to  VV 
in  Figures  3-20  and  3-21.  This  discrepancy  is  likely  due  to  a  higher  background  for  VV  compared 
to  HH.  Finally.  Figures  3-24  to  3-27  show  the  results  at  the  142°  bistatic  angle.  In  Figure  3-24, 
the  specular  RCS  HH  component  drops  by  approximately  6  dB  compared  to  the  120°  bistatic  case, 
while  the  VV  component  is  relatively  unchanged.  There  is  a  difference  of  several  decibels  in  the 
simulated  and  measured  specular  responses  in  both  HH  and  VV  components.  It  is  further  noted 
in  Figure  3-24  that  the  HH  main  lobe  has  broken  up  into  three,  while  the  VV  main  lobe  continues 
to  broaden  compared  to  the  3  =  1 2G°  case.  There  are  significant  differences  between  the  simulated 
and  measured  lobe  amplitudes  in  both  Figures  3-24  and  3-25,  but  the  general  shapes  are  in  good 
agreement.  The  differences  are  attributed  to  an  increase  in  the  background  level.  In  Figure  3-27 
the  agreement  between  the  cross-polarized  measured  data  and  simulated  data  is  good.  A  mirror 
symmetry  between  HY  and  YH  components  is  evident  here  as  should  be  the  case  for  a  symmetric 
target  [19j. 

The  specular  response  of  the  copolarized  RCS  for  the  untilted  plate  is  summarized  in  Figure  3- 
28.  The  simulated  data  indicate  that  the  specular  response  decreases  monotonically  with  increasing 
bistatic  angle  and  the  decrease  is  more  rapid  for  HH  compared  to  VV.  The  measured  data  track 
the  simulated  curves  very  well  with  the  exception  of  the  VV  component  for  0  —  142°  where  the 
measured  RCS  differs  from  the  simulation  by  about  2  dB.  To  check  the  accuracy  or  convergence  of 
the  ESP-4  simulated  RCS  patterns  at  the  3  -  142°  case,  the  rhombus  plate  was  modeled  using  a 
finer  grid  of  patches  with  a  maximum  segment  length  equal  to  0.15A.  With  this  maximum  patch 
size,  the  ESP-4  software  generated  a  17  x  16  grid  of  272  patches  for  each  polarization  with  an 
actual  patch  segment  length  equal  to  0.143A.  The  total  number  of  modes  is  544,  which  is  1.74 
times  the  number  of  modes  used  with  the  0.2A  maximum  patch  size  simulations.  A  comparison 
of  the  copolarized  RCS  for  HH  and  VV  components  with  0.187A  and  0.143A  patches  (rounded 
to  0.19A  and  0.14A,  respectively)  for  the  3  =  142°  untilted  plate  is  made  in  Figure  3-29.  The 
simulated  specular  values  agree  to  within  0.1  dB.  and  the  overall  good  agreement  between  the 
patterns  indicates  an  insensitivity  to  the  patch  size,  Thus,  the  differences  between  the  measured 
and  simulated  specular  responses  (as  well  as  the  responses  at  other  angles)  observed  in  Figure  3-28 
are  attributed  to  measurement  error. 
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ASPECT  ANGLE  9(d*fl) 


Figure  3-10.  Comparison  of  ESP-4  simulations  and  measurements  of  the  copolarized 
(HH,  VV)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic  angle  is 
0  =  15°  and  there  is  no  plate  tilt. 
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Figure  3-12,  Comparison  of  ESP-4  simulations  and  measurements  of  the  copolarized 
(HH,  VV)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistalic  angle  is 
0  ~  45°  and  there  is  no  plate  tilt. 
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Figure  3-14.  Comparison  of  ESP-4  simulations  and  measurements  of  the  cross- 
polarized  (HV,  VH)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic 
angle  is  0  -  45°,  and  there  is  no  plate  tilt. 
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Figure  3-15,  Comparison  of  ESP-4  simulations  and  measurements  of  tire  cross- 
pol  seized  (HV,  VH)  bistatic  RCS  for  the  35-  x  27-in  inch  plate  at  1,3  GHz,  The 
bistatic  angle  is  0  —  45' ,  and  the  plate  is  tilted  45°. 
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Figure  3-16,  Comparison  of  ESP-4  simulations  and  measurements  of  the  copolarized 
(HH,  VV)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic  angle  is 
0  =  90®,  and  there  is  no  plate  tilt. 
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Figure  3-1 7.  Comparison  of  ESP-4  simulations  and  measurements  of  the  copolarized 
( HH ,  VV)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic  angle  is 
0  =  90°,  and  the  plate  is  tilted  45°, 
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Figure  3-18.  Comparison  of  ESP-4  simulations  and  measurements  of  the  cross- 
polarized  (HV,  VH)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic 
angle  is  0-90°,  and  there  is  no  plate  tilt. 
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Figure  3-19.  Comparison  of  ESP-4  simulations  and  measurements  of  the  cross- 
polarized  (HV,  VH)  bistatic  RCS  for  the  35 -  x  27-in  plate  at  1.3  GHz.  The  bistatic 
angle  is  0  =  90°,  and  the  plate  is  tilted  45c‘. 
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Figure  3-20 .  Comparison  of  ESP-4  simulations  and  measurements  of  the  copolarized 
(HH,  W)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic  angle  is 
0  =  120%  and  there  is  no  plate  tilt. 
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Figure  3-21.  Comparison  of  ESP-4  simulations  and  measurements  of  the  copolarized 
(HH,  VV)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic  angle  is 
0  —  120',  and  the  plate  is  tilted  45 
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Figure  3-23.  Comparison  of  ESP-4  simulations  and  measurements  of  the  cross- 
polarized  (HV,  VH)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic 
angle  is  0  =  120°,  and  the  plate  is  tilted  45°. 
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Figure  3-24.  Comparison  of  ESP-4  simuiations  and  measurements  of  the  copolarized 
(HH,  V\r)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic  angle  is 
0  -  142°,  and  there  is  no  plate  tilt. 
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Figure  3-25.  Comparison  of  ESP-4  simulations  and  measurements  of  the  copolarized 
(HH,  W)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic  angle  is 
0  -  142°,  and  the  plate  is  tilted  45\ 
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Figure  3-26.  Comparison  of  ESP-4  simulations  and  measurements  of  the  cross- 
polarized  (HV  VH)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic 
angle  is  0  —  142° ,  and  there  is  no  plate  tilt. 


34 


0  90  180  270  360 

ASPECT  ANGLE  0(d«o) 

Figure  3-27.  Comparison  of  ESP-4  simulations  and  measurements  of  the  cross- 
polarized  (HV,  VH)  histatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz.  The  bistatic 
angle  is  0  —  142°,  and  the  plate  is  tilted  45°. 
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Figure  3-28 .  Comparison  of  ESP-4  simulations  and  measurements  of  the  specular 
copolarized  (HH,  VV)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz. 
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Figure  3-29.  Convergence  check  of  ESP-4  simulations.  The  figure  shows  the  copolar¬ 
ized  (HH,  VV)  bistatic  RCS  for  the  35-  x  27-in  plate  at  1.3  GHz  with  two  different 
surface  patch  sizes.  The  bistatic  angle  is  f3  —  142°,  and  there  is  no  plate  tilt. 
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4.  CONCLUSION 


This  report  has  described  modifications  made  at  Lincoln  Laboratory  to  The  Ohio  State  Uni¬ 
versity  electromagnetic  surface  patch  code  version  4  (ESP-4)  and  its  utilization  in  the  computation 
of  the  bistatic  radar  cross  section  of  a  rhombus-shaped  flat  plate.  The  modifications  are  the  re¬ 
structuring  of  the  input  data  statements  and  the  addition  of  a  new  option  for  bistatic  radar  cross 
section  computation  with  a  fixed  bistatic  angle,  A  listing  of  the  revised  software  is  given  in  the 
appendix. 

Comparisons  of  the  ESP-4  method  of  moments  simulations  with  far-field  measurements  over 
bistatic  angles  15°  to  142°  have  been  made  at  1.3  GHz  for  a  35-  x  27-in  rhombus  plate  in  untilted 
and  tilted  configurations  and  good  agreement  is  observed.  From  these  simulations,  it  is  concluded 
that  both  the  copolarized  and  cross-polarized  bistatic  RCS  of  arbitrarily  oriented  flat  plates  can  be 
accurately  predicted  with  this  software. 
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APPENDIX  A 

REVISED  ESP-4  SOFTWARE  LISTING 


The  purpose  of  this  appendix  is  to  list  the  main  program  of  a  modified  version  of  the  electro¬ 
magnetic  surface  patch  code  (ESP-4).  This  software  was  obtained  in  July  1988  from  The  Ohio  State 
University  and  has  since  been  modified  at  Lincoln  Laboratory.  The  important  modifications  are 
the  addition  of  namelist  input  data  and  the  option  for  bistatic  radar  cross  section  calculation  with  a 
fixed  bistatic  angle  and  variable  target  rotation,  Other  minor  changes  to  the  code  are  documented 
with  comments  in  the  modified  ESP-4  code  listed  below.  The  main  program  contains  nearly  1700 
lines  of  code  and  is  listed  in  its  entirety.  The  ESP-4  subroutines  contain  an  additional  7000  lines  of 
code  and  are  not  listed  here.  However,  in  the  subroutines  no  changes  other  than  a  global  conversion 
from  (ARCOS.  ARSIN'.  ARTAN)  to  (ACOS.  ASIN.  ATAN)  have  been  made.  Following  the  listing 
are  two  data  files  corresponding  to  the  untilted  and  tilted  plate  configurations. 
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•♦•♦•Modlf  lad  ESP-4  FORTRAN  77  cod*  liit.inj  (SUN  3  eoaputar  ayatas) . 
•  •••Hair.  progran  . ilanaaa ;  aap4nam,T  at  Lincoln  Laboratory 


C  ESP00010 
C  ESP00020 
C  OHIO  STATE  UNIVERSITY  ELECTROSCIENCE  LAB  MOMENT  METHOD,  SURFACE  ESP00030 
C  PATCH  CODE.  BASED  UPON  PIECEWISE  SINUSOIDAL  REACTION  ESP00C40 
C  FORMULATION.  TREATS:  ESP00060 
C  THIN  WIRES  ESP00060 
C  RECTANGULAR  OR  POLYGONAL  PLATES  ESP00073 
C  WIRE/PLATE  ATTACHMENTS  (AT  LEAST  0,1. WAVE  FROM  EDGE)  ESP00080 
C  PLATE/PLATE  ATTACHMENTS  ESP00090 
C  OPEN  OR  CLOSED  SURFACES  ESP00100 
C  E5P00110 


C  REFERENCE  "A  USER'S  MANUAL  FOR  AS  ELECTROMAGNETIC  SURFACE 

C  PATCH  CODE:  ESP  VERSION  IV,"  BY  E.H,  NEWMAN,  OSU/F.SL  REPORT 

C  718199-11 ,  AUGUST  1988, 

C 

C  ALSO  "A  USER’S  MANUAL  FOR  ELECTROMAGNETIC  SURFACE  PATCH  CODE:  'ESP00120 

C  VERSION  III  -  POLYGONAL  PLATES  AND  WIRES"  BY  E.H.  NEWMAN  AND  E3P00130 

C  R.L.  DILSAVO'  OSU/ESL  REPORT  716148-19,  APRIL  1987.  ESP00140 

C  ESPOO 1  BO 


C  ANY  COMMENTS  CAN  BE  REFERRED  TO: 


ESPOO 160 


C 

C 

C 

C 

C 

C 

C 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 


ESP00170 

EDWARD  H.  NEWMAN 

ESPOOlBO 

ELECTROSCIENCE  LABORATORY 

ESP00190 

1320  KINNEAR  ROAD 

ESP00200 

COLUMBUS,  OHIO  43212 

ESP’00210 

PHONE:  (614)  292-4999 

ESP00220 

ESP00230 

FILENAME  aap4naiu.f  SUN  3  COMPUTER  SYSTEM  at  MIT  LINCOLN  LABORATORY. 

MODIFIED  FOR  NAMELIST  INPUT  DATA  AND  FIXED 
BISTATIC  ANGLE  WITH  VARIABLE  TARGET  ROTATION. 

ALL  LINCOLN  LABORATORY  MODIFICATIONS  ARE 
MADE  EVIDENT  BY  THE  ABSENCE  OF  ESP  LINE  NUMBERS. 
THE  SUBROUTINES  USED  BY  aap4naa.l  HAVE  BEEN 
DIVIDED  INTO  TWO  FILES  a«p4aubal,l  AND 
«>P4aub(2.^  WHICH  ARE  THE  ORIGINAL  ESP4 
SUBROUTINES  (I.E.  NO  MODIFICATIONS). 

DATE:  3  JANUARY  1090 


C  ESP00260 
C  ESP00260 
C  ESP00270 
C  DIMENSION  INDICATORS:  ESPOO2B0 
C  THE  FOLLOWING  ARE  DEFINED  BY  THE  USER:  ESP00290 
C  ESP00300 
C  IDFIL;  INDICATOR  TO  DIMENSION  FOR  FILAMENT  TESTING  ESP00310 
C  ■  1  IMPLIES  DIMENSION  FOR  FILAMENT  TESTING  ESP00320 
C  >0  IMPLIES  DO  NOT  DIMENSION  FOR  FILAMENT  TESTING  ESP00330 
C  1DSUR:  INDICATOR  TO  DIMENSION  FOR  FULL  SURFACE  TESTING  E5F00340 
C  -  1  IMPLIES  DIMENSION  FOR  FULL  SURFACE  TESTING  ESP003S0 
C  >0  IMPLIES  DO  NOT  DIMENSION  FOR  FULL  SURFACE  TESTING  ESP00360 
C  IDWR-  MAX.  NUMBER  OF  WIRE  POINTS,  WIRE  SEGMENTS,  AND  WIRE  MODES  ESP00370 
C  IPL“  MAX.  NUMBER  PLATES  ESP00380 
C  ICN-  MAX.  NUMBER  OF  CORNERS  ON  POLYGONAL  PLATES  ESP00390 
C  IAT«  MAX.  NUMBER  WIRE/PLATE  ATTACHMENT  POINTS  ESP00400 
C  I TOT*  MAX.  TOTAL  NUMBER  OF  MODES  (WIRE*PLATE+ATTCH . )  ESP00410 
C  1DMZI  *  INDICATOR  FOR  MAKING  FREQUENCY  SWEEP  COMPUTATION  E3F00420 
C  *  1  IMPLIES  SET  UP  Z  ARRAY  FOR  INTERPOLATION  ESP00430 
0  *  0  NO  FREQUENCY  SWEEP  ESP00440 


C  IERVSR  *  MAX.  I  FOR  THE  ERVSR(IAT.I)  ARRAY  USED  IN  SUB.  ZTOTZ  ESP004 .0 
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c 

ESP00460 

c 

TME  FOLLOWING  ARE  COMPUTED  BY  THE  CODEi 

ESP00470 

c 

ISP004SO 

c 

IDZT**  MAX.  LENGTH  OF  1-D  IMPEDANCE  ARRAY  2T 

ESP00490 

c 

IDZTf "MAX .  INDICATOR  FOR  2-0  ARRAY  ZTF  USED  FOR  FILAMENT  TESTINGESP00500 

c 

IDWR2*  2*1DWR 

ESP00610 

c 

ITW2-  THE  LARGER  OF  IDWR2  AND  ITOT 

ESP00620 

c 

ESP00630 

c 

DEFINE  DIMENSION  INDICATORS 

ESP0064O 

c 

ESP00B50 

PARAMETER  (IDFIL-1) 

ESP00B60 

PARAMETER  (IDSUR-1) 

ESP00B70 

PARAMETER  (IDWR-30) 

ESPOOBBO 

PARAMETER  (IPL«20) 

E3P00690 

PARAMETER  (ICN-6) 

ESP00600 

PARAMETER  UAT-4) 

ESP00610 

c 

PARAMETER  (1T0T-200) 

ESP00620 

PARAMETER  UT0T-500) 

PARAMETER  (IDMZI  ■  1) 

ESP00630 

PARAMETER  (IERVSR  -  200) 

ESP00640 

C....THE  FOLLOWING  LINE  ADDED  FOR  NAMELIST  MODIFICATION 

PARAMETER  (INFPT-SO) 

C 

ESP00650 

C 

THE  FOLLOWING  IS  EQUIVALENT  TO: 

ESP00660 

C 

PARAMETER  (IDZT»MAXO<  <I0VR**2*IDWR) /2 , IDSUR* (IT0T»»2+IT0T)/2 , 1 >  >ESP00670 

C 

ESP00680 

PARAMETER  (NUMl«l+(IDVR**2-*IDWR)/2) 

ESP00690 

PARAMETER  <mJM2«l*  IDSUR*  ( IT0T**2+IT0T) /2) 

ESP00700 

PARAMETER  (MLT1 •  ( UNUM1 /MIM2 ) *NUM2 ) / (1 . 0*NUM 1 ) ) +0 . 99099 > 

ESP007 10 

PARAMETER  (MLT2-  ( <  (NUM2/NUM1 )  *NUM1 )  /  (1 . 0*NUH2) )  +0 .  99999) 

ESP00720 

PARAMETER  (NUM3-I-1+ (MLT1*NUN1+HLT2*NUM2) / (MLT1+MLT2) ) 

ESP00730 

PARAMETER  (NUM4»1+1) 

ESP00740 

PARAMETER  (MLT3* ( ( (NUM3/NUM4 ) *NUM4 )  /  ( 1 . 0*NUM3) ) +0 . 99999 ) 

ESP007BO 

PARAMETER  (MLT4-< ( (NUM4/NUM3) *NUM3) /Cl . 0»NUH4) )»0 . 99999) 

ESP00760 

PARAMETER  (IDZT— 1-*<MLT3»NUM3*-MLT4»NUM4)/<MLT3*MLT4) ) 

ESP00770 

C 

ESP00780 

c 

THE  FOLLOWING  IS  EQUIVALENT  TO: 

ESP00790 

c 

PARAMETER  (IDZTF-MAXOUDFIL‘ITOT.1)) 

ESP00800 

c 

ESP00810 

PARAMETER  (NUM6"1+IDFIL*IT0T) 

ESPOO02O 

PARAMETER  (NUM6M+1) 

ESP00830 

PARAMETER  (KLT6-C ( (NUM6/NUM6) *NUM6)/ (1 . 0*NUM6) )+0 . 99999) 

ESP00840 

PARAMETER  (MLT6-( ( (NUM8/NUM6) »NUMS)/ < 1 . 0*NUM6) )+0 . 99999) 

ESP008EO 

PARAMETER  ( IDZTF—  1+ (MLT6*NUM5+MLT6*NUM6) / (MLT6+MLT6) ) 

ESP00860 

c 

ESP00870 

PARAMETER  (IDUR2-2*IDWR) 

ESP00880 

c 

ESP00890 

c 

THE  FOLLOWING  IS  EQUIVALENT  TO: 

KSP00900 

c 

PARAMETER  (ITV2-MAXO!ITOT,IDWR2)) 

ESP00910 

c 

ESP00920 

PARAMETER  (NUM7-1+IT0T) 

ESP00930 

PARAMETER  (NUM8-1»IDWR2) 

ESP00940 

PARAMETER  (KLT7- ( ( (NUHT/NUMB) »NUM8 ) / (1 . 0»NUM7 ) ) +0 . 99999 ) 

ESP009B0 

PARAMETER  (MLT8» ( ( (NUM3/HUM7 ) *NUM7)/ ( i . 0»NUM8) ) +0 . 99999) 

ESPOOBBO 

PARAMETER  UTW2— 1+<MLT7*NUM7+KLT8*NUM8)/ (HLT7+MLT6) > 

ESP00970 

c 

ESP00980 

c 

IDMI  ■  1  IF  IDMZI  »  0 

ESP00990 

c 

■  3  IF  IDMZI  «  1 

ESP01000 

PARAMETER  UDHI«1+IDMZI*2) 

ESP01010 

PARAMETER  <IDZTI-1+IDMZI»IDZT) 

ESP01020 
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c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PARAMETER  (IDZTFI-UIDMZI. IDZTF)  ESP01030 

ESP01040 

ESP01060 

THE  FOLLOWING  ARE  DIMENSIONED  BY  IDWR:  ESP01060 

ESP01070 

COMPLEICGD(IDUR) , SGD(IDUR)  ESPOIOBO 

DIMENSIONS  IDWR) , lA(IDVR)  ,IB<IDWR)  .NDdnWR) .  ISC(IDVR)  ,MD(IDVR,4)ESPOi090 

DIHENSIONll (IDWR)  ,I2(IDWR)  , I3CIDVR) , JA(IDWR) , JB(IDWR)  ESPOUOO 

DIMENSIONS  IDWR)  ,Y(IDWR)  ,Z(IDWP.)  ESPOlliO 

2SP01120 

THE  FOLLOWING  ARE  DIMENSIONED  BY  IDWR2  •  2*IDVR:  ESP01I30 

ESP01140 

CGMPLEXZLD ( IDWR2) ,VG(IDWR2)  ESPOIUO 

ESP01160 


THE  FOLLOWING  ARE  DIMENSIONED  BY  THE  LARGER  OF  ITOT  OR  I DVRS  ESP01170 


ESPoneo 

C0MPLEXCQ(ITW2)  ESP01190 

ESP01200 

THE  FOLLOWING  ARE  DIMENSIONED  BY  IPL  AND  ICN:  ESP01210 

ESP01220 

DIMENSI0HNM12N(IPL) ,NM23N(IPL) ,IPN(IPL) , PCN(3, ICN.IPL)  ESP01230 

DIMENSION  XP(ICN) ,YP(1CN) ,ZP(ICN) 

DIMENS10NNCNRS (IPL) , SEGMCIPL) ,PC(3,ICN) .NPLll(IPL) ,NDNPLT(IPL)  ESP01240 
DIMENSION  NPL22C2PL) ,IREC(IPL) ,SUBI(IPL) .IQSCIPL)  ESP01260 

COMPLEX  ZSHT(IPL),ZSHTF<IPL)  ESP01260 

ESP01270 

THE  FOLLOWING  ARE  DIMENSIONED  BY  I AT  AND  IERVSR:  ESP01230 

ESPOUOO 

INTEGER  NASAT(IAT) ,IABATCIAT) 

DIMENSIONBDSKCIAT) .RMIN(IAT) , DR(IAT) .NPLA(IAT) ,NSA(IAT) .  ESP01300 

2  PDIST(IAT)  ESP01310 

COMPLEXERVSR C I AT , IERVSR) .ZLDA(IAT) ,VGA(IAT)  ESP01320 

ESP01330 

THE  FOLLOWING  ARE  DIMENSIONED  BY  ITOT :  ESP01340 

ESP013S0 

DIMENS10NPA(ITOT,4i3)  .PBaTOT.4,3)  .IQUAD(ITOT)  ESP01360 

DIMENSIONIOVT(!VOT,4) ,DOVL(ITOT) , ITK(ITQT) , OVEPCITOT, 3 , 2)  ESPOI370 

COMPLEXC J ( I TOT ) .CJP(ITOT) .CJT(ITOT) .ETT(ITOT) .EPP(ITOT) ,V(ITOT)  ESPO13B0 
COMPLEX  EXN(ITOT)  .EYN(ITOT)  .EZN(ITOT)  ESP01390 

DIMENSION  MPLA(ITOT)  .MPLB(ITOT)  .PSZ(ITOT)  ESPOUOO 

ESP01410 


THE  FOLLOWING  ARE  DIMENSIONED  BY  INFPT  ••••NAMELIST  DATA  — 
INTEGER  IFMM( INFPT) .IABB(INFPT) 

COMPLEX  VLCG(INFPT) .ZLL(INFPT) 


THE  FOLLOWING  AAE  DIMENSIONED  BY  IDZT:  ESP01420 

ESP01430 

COMPLEXZT (IDZT)  ESP01440 

ESP01460 

THE  FOLLOWING  ARE  DIMENSIONED  BY  IDZTF;  ESP01460 

ESPO1470 

COMPLEX  ZTFUDZTF,  IDZTF)  ESP014B0 

espoubo 

THE  FOLLOWING  ARRAYS  HOLD  THE  Z  MATRICES  FOR  THE  FREQ.  SWEEP  ESPOUOO 

ESPOUIO 

COMPLEX  ZTIN(IDMI,IDZTI) ,ZTFTN(IDMI .IDZTFI ,IDZTFI)  ESP01S20 

ESP01630 

THE  FOLLOWING  ARE  FIXED  DIMENSIONED:  ESP01640 

ESP01S60 

DIMENSION  2ZC(3)  ESPOUBO 
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nei  n  n  anon  non 


C0MPLEXETEU441)  ,£P£(1441)  ,ERE(144l)  ,£142(1441)  , 

2  EPAZU441) 

DHffiNSION  PET(1441) ,PEP(1441) ,PER(1441) 

COMPLEX  ERAZU441) 

COMPLEXET* , 0AM , 2 1 1 , Y 1 1 , ERRS , ETTS , EPPS , ETPS , EPtS , ZS , ZL , VLQ , ECD 
COMPLEX  ZIN.YIN.VIN 

COMPLEXCOA , CGB.CGCN , CEXP 1 , CEXPS , CEXP3 , CEXPK ,FACl ,FAC2 , FACS , PACK 
COMPLEX  EXT , EYt , EZT , EXP , EYP , E2P , ETRS , EPRS , X  J , YSHT , YSHTF 
COMMON  /A/  WV,PI,A.Q,OAM,ETA,XK 

c»**add:tion  of  namelists 

NAMELIST  /RNCm/NOO ,  NPRINT ,  NRUNS ,  NWOS ,  I«R ,  I WRZT ,  I  NT ,  1NTP ,  INTO , 
2INWR , IRGM i IFIL ,RF , INDZI 

MAMELIST  /FSWEEP/FMC1 , FMC2 , DF2I , OFF , IRS 12 , THRD , PHRD , THRI , PHR1 
NAMELIST  /PATTRN/IFE , IPFE , FNDFE , PHFE , IFA , IPFA ,FNDFA , THFA , 

2ISE , IPSE , FNDSE ,PKSE ,THIN ,PHIN , ISA , IPSA , FNDSA , THSA , 

3AZRANQ .ELRANG , AZMIN , ELMIN , IBISC , BETA , NPTBIS , BANGRG 
C--NOTE:  AZRANG ,  EL RANG  INCLUDED  TO  REPLACE  FIXED  360  DEC. 

C*»NOTE:  AZMIN, ELMIN  INCLUDED  TO  START  PATTERNS  AT  ARB.  POINT 
C--NOTE:  18180*1  USES  FIXED  BISTATIC  ANGLE  FEATURE  WITH 
C  BETA  THE  BISTATIC  ANGLE, 

C  NPTBIS  THE  NUMBER  OF  TARGET  ROTATION  ANGLES, 

C  BANGRG  THE  ANGULAR  RANGE  OF  TARCET  ROTATION 

C  NOTE:  BANGRG  CAN  BE  >0  OR  <0  TO  DETERMINE  ROTATION  DIRECTION 

C--WHEN  THE  BISTATIC  OPTION  IS  USED,  ELRANG  AND  AZRANG  ARE  SET-0 
NAMELIST  /FVIRET/FNC , CMM , A , NPLTS 

NAMELIST  /PLATEG/NCNRS , SEGM , IREC . I PN , 3 GS , ZSHT , XP , YP , ZP 
NAMELIST  /SAVEZ/IWRZM.IRDZM 

NAMELIST  /WIREAG/NM ,NP , NAT , NFPT , NFS  1 , NFS2 ,X,Y,Z,IA,IB 
NAMELIST  /GENLOD/ I FMM , I ABB , VLGG , ZLL 
NAMELIST  /ATTACH/NASAT . IABAT , NPLA , VGA , ZLDA , BOSK 
CHARACTER  DATNAM*25,0UTNAM»26 
•  ' J(I ,J ,NTOT)-( J-1)-NT0T- (J»J-J)/2*I 
AMP(GAM) -CABS (GAM) 

C  PHS(GAM)-180 ,0»BTAN2(AIMAG(0AM) , REAL (GAM) )/3. 141B926 

C  OPEN (UNIT-6,  NAME- 1 OUTFL . DAT 1 , TYPE- ' UNKNOWN ’ ) 

C 

C  THE  MAIN  PROGRAM  CONTAINS  FOUR  CALLS  TO  THE  CLOCK  FUNCTION 

C  CETCP(I),  WHERE  I  IS  THE  CLOCK  READING  IN  HUNDREDTHS  OF  A 

C  SECOND.  THESE  FOUR  LINES  HAVE  BEEN  "COMMENTED"  OUT.  IN  ORDER 

C  FOR  THE  CODE  TO  OUTPUT  CPU  TIMES,  THE  USER  MUST  REPLACE  THESE 

CALLS  TO  GETCP  BY  A  COMPARABLE  CLOCK  FUNCTION  ON  HIS  SYSTEM. 

CALL  GETCP2  PERFORMS  THE  CPU  TIME  FUNCTION  ON  THE  SUN  3  SYSTEM 
CALL  GETCP2 (ICPU) 

PI-3 . 141 69265 
19999-0 
XJ-(0. 0,1.0) 

IANT-0 

INM-  MAX.  NUMBER  WIRE  SEGMENTS 
ICJ-  MAX.  NUMBER  WIRE  MODES 
XPLM-  MAX .  NUMBER  PLATE  MODES 
ICC-  MAX.  SIZE  OF  2-0  ARRAY  ZTF 
IHM-IDWR 
ICJ-IDWR 
IPLM-ITOT 
ICC-IDZTF 

NGO  -  0  IMPLIES  SET  UP  GEOMETRY  BUT  DO  NOT  RUN. 

NGO  -  1  IMPLIES  RUN. 


ESP0167C 

ESF01660 

ESP01E90 

ESP0160C 

ESP01610 

ESP01620 

ESP01630 

E3P01640 

ESP01660 


ESP01660 

ESP01670 

ESP01680 

ESP016W0 

ESP01700 

ESP01710 

ESP01720 

ESPO1730 

ESP01740 

ESP01760 

ESP01760 

ESP01770 

ESP01780 

ESP01790 

ESP01800 

Espoieio 

ESP01820 

ESP01830 

ESP01840 

ESP018S0 

ESP01860 

ESP01870 

ESP01880 

ESP01890 

ESP01900 

ESP01910 

ESP01920 

ESP01930 
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C  NPRINT  •  0  IMPLIES  PRINT  INPUT  DATA  ONLY  ESP01940 

C  NPRIHT  »  1  IMPLIES  PRINT  WIRE  AND  PLATE  GEOMETRY .  ESP019B0 

C  MPRINT  »  2  IMPLIES  PRINT  INPUT  DATA  AND  WIRE/PLATE  GEOMETRY.  ESP01960 

C  NPRINT  ■  3  IMPLIES  PRINT  NEITHER.  ISP019T0 

C  ESP01980 

WR1TEC6, 29691 


2969  FORMAT  OX, 'ENTER  INPUT  DATA  PILE  NAME  (typically  **p4ata,datajO  1 ) 
READ(6,«)DATNAM 

OPEN ( 1 1 . FI LE-DATN AM , FORM* > FORMATTED  > ) 

URITEC6, 39691 

3969  FORMATdX, 'ENTER  OUTPUT  DATA  FILE  NAME  (typ.  «ip4nasf ort . 6)  >) 
READ(5,»)0UTNAM 

OPEN (8 , FILE-OUTNAM , FORM- ' FORMATTED ' ) 

C  READtll, *1X00, NPRINT, NRUNS.NWGS.IWn.IWRZT.INT.INTP, INTO, INWR.IRGM.ESP01990 


C  1IF1L.RF, INDZI  ESP02000 

READ (11,RN CTRL) 

C  ESP02010 

C  READ  IN  PARAMETERS  OF  FREQUENCY  SWEEP  COMPUTATION  ESP02020 

C  ESP02030 

12468-0  ESP02040 

IFCINDZI . NE.01TKEN  ESP020S0 

IFdDMZl  .EQ.01TKEN  ESP02000 

VRITE(6,3216)  ESP020T0 

3216  F0RMAK3X,  'DIMENSION  INDICATOR  IDMZI  MUST  BE  SET  TO  1  IF  INDZI  '/ESP02080 

2  31 . 'IN  READ  1  IS  SET  TO  1  OR  2'}  ESP02090 

STOP  ESP02100 

ENDIF  ESP02110 

C  READdl  ,*)FMC1  ,FMC2,DFZI  ,DFF,  IRS12  ,THRD,PHRD,THINC,PHINC  ESP02120 

READ(ll.FSWEEP) 
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noon 


DPTOT- (FMC2-FMG 1 ) 

MPZI »0 . 99+DPT0T/DFZI 
IF(NFZI.LT.2)NFZI«2 
BFZI-DFTOT/NFZI 
NFZI-NF22+1 
NFF-0,99*CFT0T/DFF 
IFCMFF.LT, i)NFF*l 
DFF-DFTOT/NFF 
MFF-NFF+ 1 

:FCNOO.ME,0)WRITE(10,3210)RF,THRL‘>PHRL,THIMC,PHHIC11RS12, 

2  MFFtNFZI , INDZI 

3210  FQRMAT(1X,E11,6,4(1X,F6.1),1X,6I4) 

II2LST-2 
FNDFE- 1 . 0 
FNDFA-1,0 
ENDS';- 1.0 
FNDSA-1.0 
END  IF 

IF(INVR.EQ.O)THEN 

NM“0 

NP-0 

NAT-0 

NVR-0 

NFPT-0 

NFS1-0 

MFS2-0 

END  IF 

IP(INT.GTiO)INT-2*((INT+l)/2) 

JNTP-2* ( ( INTF+ 1 ) /2) 

INTD-2*  ( ClNTD+D/2) 

270  FGRMATC2I1 ,213) 

C 

C  READ  IN  PARAMETERS  OF  ELEVATION  AND  AZIMUTH  PATTERNS 

C  SET  ALL  TO  0  IF  INDZI  >  0 

C 

IF (INDZI . EQ . 0) THEN 
READ  ( 1 1 ,  • )  IFE .  IPFE ,  FNDFE ,  PHFE 
READ ( 1 1 , • ) IFA , IPFA . FNDFA , THFA 
READ ( 1 1 . • ) I3E , IPSE . FNDSE . PHSE . THIN . PHIN 
READ ( 1 1 ,*) ISA , IPSA . FNDSA , THSA 
CC*«*DEFAULT  NAMELIST  FATTEN  VALUES  FOR  AZRANC .ELRANO, AZMIN ,ELMIN , . . , 
AZRANG-360 . 

ELRANG-360 . 

AZMIN-0 . 

ELMIN-O. 


ESP02130 

ESP02140 

ESP02160 

ESP02160 

ESP02170 

ESP02100 

ESP02190 

ESP02200 

ESP02210 

ESP02220 

EDP02230 

E3P02240 

ESP022S0 

ESP02260 

ESPO2270 

ESP02280 

ESP01290 

ESP02300 

ESP02310 

ESP02320 

ESP02330 

ESP02340 

ESP023S0 

ESP02360 

ESP02370 

ESP02380 

ESP02390 

ESP02400 

E3P02410 

ESP02420 

ESP02430 

ESP02440 

ESP02460 

ESP02460 

ESP02470 

ESP02480 

ESP02490 

ESP02S00 

ESP02S10 

ESP02B20 


IBISC-0 

BETA-0. 

MPTBIS-1 

BANCRG-O. 

READ (11, PATTRN ) 
IFCIBISC.EQ .0)NFTBIS«1 
IFCIBISC.EQ .  DELRANG-0 . 


ELSE 

ESP02S30 

IFE-0 

ESP02S40 

IPFE-0 

ESP02660 

FNDFE- 1.0 

ESP02S60 

PHFEjO-O 

ESP02E70 

IFA-0 

ESP02680 

IPFA-0 

ESP02S90 

FNDFA-1 .0 

ESP02600 
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«.*  a 


THFA-0,0 

ISE“0 

XPSE-0 

FND3E-1.0 

PMSE-0,0 

THIN-0,0 

PH1N-0.0 

ISA-0 

IP8A-0 

FNDSA-1.0 

THSA-O.O 

ENDIF 

312  F0RMATC3I2 ,F10 . 8) 

XMAOE-O 

IFdND7I.EQ,0)THEN 

IF ( !S£ . LT . 0 , OR . ISA . LT . 0) IMAOE-i 

ELSE 

IFdRS12,LT,0)IMAGE-l 

ENDIF 

ISE*IABS(1SE) 

ISA-IABSdSA) 

IFF-0 

XSCAT-0 

I'  INDZI.EQ.O)THEN 

IFF*-]  FA+ISE+ISA  .(IE.  l)  IFF-l 
IFdSE.  EQ ,  1 .  UR ,  ISA .  EQ .  1 )  I  SCAT- 1 
IF dSE . EQ . 2 . OR , ISA . EQ . 2) ISCAT-2 
IF (ISE . EQ . 3 . OR , ISA . EQ . 3) ISCAT-3 
ELSE 
IFF-1 

IF(IABS(IRS1I) ,U].2)THEN 
ISCAT-2 

IF( ABS(THINC-THRD) *ABS (PHINC-PHRD ) . LT . 0 . 001) ISCAT-1 

ISE-1 

ISA-0 

ENDIF 

IF(IABSdRS12)  .EQ.  DTMEN 

I3CAT-0 

IFE-1 

IFA-0 

ENDIF 

ENDIF 

C  NPLOTS  -  THE  NUMBER  OF  PATTERN  PLOTS 

C  IRS12  -  1  OR  2  FOR  RADIATION  OR  SCATTERING  PATTERNS 

MPLtlTS-0 

IF (ISCAT . EQ . 0) THEN 

IF  UFE . NE . 0 . AND . IPFE . EQ . 1 > NPLOTS-NPLOTS + 1 
IF  dFA .  NE .  0 .  AND .  IPFA .  EQ .  1 )  NPLOTS-NPLOTS+ 1 
ELSE 

IF (ISE . NE . 0 . AND . IPSE . EQ . 1 ) NPLOTS-NPLOTS+ 1 
IFdSA .  NE .  0 .  AND .  IPSA .  EQ .  1 )  NPLOTS-NPLOTS*  1 
ENDIF 

IFdttDZI  .EQ.O)TH£N 
IRS12-1 

IF(I3CAT.NE.0)IRS12-2 

ENDIF 

DO  700  NRUN-l.NRUNS 

READ  FREQUENCY,  WIRE  CONDUCTIVITY,  AND  WIRE  RADIUS 


ESP02610 

ESP02620 

ESP02630 

ESP02640 

ESP02660 

ESPO2660 

ESP02670 

ESP02680 

ESP02690 

ESP02700 

ESP02710 

ESP02720 

ESP02730 

ESP02740 

ESP02760 

ESP02760 

ESP02770 

ESP02780 

ESP027P0 

ESP02800 

ESP02810 

ESP02820 

ESP02830 

ESP02840 

ESPOlABO 

ESP02800 

ESP02870 

ESP02880 

ESP02890 

ESP02900 

ESP02910 

ESP02B20 

ESPG2S30 

ESPC2940 

ESP029BD 

ESP02960 

ESP0297C 

ESP029B0 

ESP02990 

ESP03000 

ESP03010 

ESP03020 

ESP03030 

ESP03040 

ESP03060 

ESP03060 

ESP03070 

ESP03080 

ESP03090 

E8P03100 

ESP03110 

ESP03120 

ESP0U130 

ESP03140 

ESP031SO 

ESP03160 

ESPO:i70 

ESP03180 

ESP031W0 

ESP03200 
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C  READCll, OFHC.CMM, A 

READCll, FWIRET) 

IF  < 1NDZI , NE . 0)FMC*KMC2 

IFCNPLOTS . OT . 0 , AND . INDZt .  EQ . 6 . AHD . NOO . HE , 0) 

2  WRITE(B,33?)NPL0TS, IRS12 ,FMC,RF 
337  F0RMATC1X ,2(13, IX) ,2(F11 ,4 ,2X) ) 

W-300.0/FMC 

TOUCH»O.OOt»WV 

IOK-1 

c 

C  READ  IN  THE  NUMBER  OF  PLATES 

C 

C  READCll,*) NPLTS 

6E0HX— 1.0 
IF (NPLTS .GT. IPL) THEN 
WRITE (6 ,330) NPLTS 

330  FORMAT ( 1  *•••••  INCREASE  PARAMETER  IPL  TO  *,I4) 

STOP 

END  IF 

NOVTeO 

NPLTM-0 

IF(NPLTS .  EQ. O)G0T0462 

IOKT-1 

C 

C  IF  NPLTS  <  0,  CEHERATE  PLATE  GEOMETRY  IN  COEOM 

C 

IFCNPLTS.LT. OTHEN 
C 

C  PLATE  GEOMETRY  GENERATED  IN  SUBROUTINE  CALL  BELOW 

C  SUBROUTINE  PQEOM  IS  CONTAINED  IN  ESP4PGH  FORTRAN . 

C 

CALL  PGEOM ( IPL . ICN , NPLTS , NCNRS , SEOM , IREC . IPN . I5S , ZSHT ,PCN) 

IF  (NPLTS.  GT.  IPL)  THEN 
WRITE (6 , 330) NPLTS 
STOP 
EJIDIF 
ELSE 
C 

C  READ  IN  THE  PLATE  GEOMETRY 

C 

DO  464  NPL-1, NPLTS 

C  READCll , «)NCNRS(NPL) .SEOMCNPL) .XREC(NPL) .IPN(NPL) ,IOS(NPL> , 

C  2  ZSHT(NPL) 

READCll, PLATEG) 

WRITE(e.PLATEG) 

IF ( SEGMCNPL) . GT . SEGMX) SEGMX»3E3M ( NPL ) 

IFCNCNRSCNPL) .EQ.4)I0S(NFL)>0 
IFCNCNRS(NPL) .GT.ICN)THEN 
write (6 , 331)ncnrs(npl) 

331  FORMAT ( '  INCREASE  PARAMETER  ICN  TO  AT  LEAST  •,!<) 

STOP 

ENDIF 

DO  466  NCNR»1,NCNRS(NPL) 

C  READ(11,»)PCNC1,NCNR|HPL) ,PCN(2,NCNR,NPL)  ,PCN(3,NCNR,NPL.' 

C*«*PIU-IN  PCN  ARRAY  FROM  NAKELIST  DATA 
PUN ( 1 , NC«R , NPL) «XP (NCNR ) 

PCN (2 ,NCNR , NPL) *YP (NCNR) 

PCN(3,NCNR,NPL)-ZP<NCNR) 

466  CONTINUE 
464  CONTINUE 
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END  IF 

DU  466  NPL- 1 .NPLTS 

CALL  PLPLCKtPCN , IPL , ICN , NCNRS(NPL) , TOUCH , NPL , I OK ) 

IF  (lOK.Efl.O)  S0XT“0 
465  CONTINUE 

IF(IOKT.EQ.O)  GOTO  9374 
DO  46?  NPL-1  .NPLTS 
NCNS-NCNRS(NPL) 

SEO-SEQM(NPL) 

IRE-IREC(NPL) 

IP-IPN(NPL) 

lO-IOS(NPL) 

DO  466  NC-l.NCNS 
DO  468  1-1,3 

468  PC(I,NC)-PCN(I,NC(NPL) 

C 

C  SEGMENT  PLATES  INTO  MODES 

C 

CALL  PLATE3 <  PC , NCNS , ICN , NPL , NDNPLT , PA , PB , 1PLM , SEO . 

1  1QUAD , WV , IRE , IP , MPL1 , MPL2 ,  IOK , NM12 , NM23 , IG) 

1F(NPL .  Et) .  NPLTS )  NPLTM-NDNPLT (NPLTS ) 

NPLll(NPL)-MPLl 
NPL22 (NPL) -MPL2 


NH12N(NPL)-NM12 
NM23N (NPL) -NM23 

46?  CONTINUE 
C 

C  CONNECT  TOUCHING  PLATES  WITH  OVERLAP  MODES 

C 

CALL  POPLOV ( NPLTS , PCN , NCNRS , TOUCH , 

A  SEGM , PA ,PB , NOVT , NPLTM , IPL , IPLM , ICN , IOVT , DOVL , ITK , NOPL . 
ft  TQUAD.WV, NDNPLT, OVEP) 

462  CONTINUE  . 

NPLTM-NPLTM+NOVT 
DO  600  NWC-l.NWCS 
112-1 

CALL  GETCP2(JCPU) 

IWRZM-0 

IRDZM-0 

C 

C  READ  IN  Z  MATRIX  READ/WRITE  FROM  DISK  INDICATORS 

C 

C  READ ( 1 1 , » ) I VR2M , IRDZM 

READ (11 .SAVEZ) 

C  IF'INDZI .GT.O)THEN 

C  IWRZM-0 

C  IRDZM-0 

C  ENDIF 

W-300.0/FMC 
INM2-2-INM 
DO  2774  I-1.INH2 
ZLD(1)-CMPLX<0. 0,0.0) 

VG(I)-CMPLX(0. 0,0.0) 

IF(I .CT. IAT)G0TO2?74 
V«A(I)-CMPLX(0. 0,0.0) 

ZLDA(I)-CNPU<(0. 0,0,0) 

2774  CONTINUE 

IFdNVR.EQ.O)  NWR-0 
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tmmm.iQ.o)  nat-o 

IF ( 1NWR . tQ . 0) G0T02773 
IFdROM.EQ. 0)00102800 
C 

C  SET  UP  POINTS  AND  SEGMENTS. 

C 

C  READ ( 1 1 , • ) NM , NP , N AT , NFPT , NFS1 , NFS2 

READdl.VIREAG) 

IF(NAT.OT. IAT)THEN 
WR1TE(6,332)NAT 

332  FORMAT ( ’  INCREASE  PARAMETER  I AT  TO  AT  LEAST  ‘,14 ) 

STOP 
END  IF 

C 

C  READ  IN  COOR,  OF  NP  POINTS 

C 

DO  3810  I-l.NP 

C  READ(U,»)X<I)  ,Y(I)  ,2(1) 

2810  CONTINUE 
C 

C  READ  IN  ENDPOINTS  OF  NM  SEGMENTS 

C 

DO  2820  I-l.NK 

C  READ(11 , <)IA (I) i 16(1) 

2820  CONTINUE 
C 

C  READ  IN  WIRE  GENERATORS  AND  LOADS 

C 

READdl  .CENLOD) 

DO  2830  I-l.NFPT 

C  IF(NFPT.QE.l)READdl,*)IFM,IAB1VLG,2L 

C  II-IFM+IAB»NM 

II-IFMM(I)+IABB(I).NM 

VG«I)-VLGO(I) 

2LDCII)“ZLL(I) 

2830  CONTINUE 

C 

C  READ  IN  ATTACHMENT  POINT  GEOMETRY 

C 

IF ( HAT . EQ . 0 ) G0T028S0 
READdl,  ATTACH) 

DO  2840  I-l.NAT 

C  READdl  i ONAS,  1AB ,NPLA(I)  .VGA Cl)  ,ZLDA  d)  ,BDSK(I) 

C  NSA(I)-NAS+IAB*NM 

NSA ( I ) “NASA! (I)tlABAT(I) *NM 
2840  CONTINUE 

GOT02860 
C 

C  IF  IRCH  -  0,  GENERATE  WIRE  GEOMETRY  IN  SUBROUTINE  WOEOM 

C 

2800  CALL  UGEOMdA.IB.X.Y.'Z.NM.NP.NAT.NSA.NPLA.VOA.BDSK, 

2  ZLDA.NWG, VG,2LD|WV,NFS1 ,NFS2) 

2860  CONTINUE 

IF  UFIL .  EQ .  1 .  AND ,  NAT .  GT .  0)  THEN 
WRITE(«,2862) 

2862  FORMAT (/3S, 'ERROR!  FILAMENT  TESTING  CAN  NOT  BE  USED  IF  >, 

2  ‘GEOMETRY  CONTAINS  A  WIRE/PLATE  ATTACHMENT ’ /3X , 

3  ‘SET  IFIL  •  0  IN  READ  1') 

STOP 

ENDIF 
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PLACE  HOPES  ON  THE  WINE 

CALL  SORT (IA,IB,Il,I2,I3,Jt,JB,MD,ND,NM,NP, HWR , HAX , HIN , IC J , INM) 
NDWR«MAXO(NM,NP,NWR) 

IF (NDUR i CT . IDWR) THEN 
WRITE (0, 336) NDWR 

336  FORMAT ( 1  ••••••  INCREASE  PARAMETER  IDWR  TO  AT  LEAST  ',14) 

STOP 
ENDIF 
CONTINUE 


2773 


DETERMINE  ON  WHICH  PLATES  SURFACE  PATCH  MODES  LIE 


DO  2048  I-J.NPLTM 

CALL  WCHPLT ( I , NPLTS , NOPL , NDNPLT , 1 0  VT , ITK . ITOT , MPLA , MPLB ) 
204S  CONTINUE 

NTOT-NWR+NPLTM+NAT 

IF(NTOT.OT.ITOT)THEN 

WR1TE(6,334)NT0T 

334  FORMAT ( 1  ••*•••  INCREASE  PARAMETER  ITOT  TO  AT  LEAST  • .IB) 
STOP 
ENDIF 

N1«NPLTM-N0VT+1 

N2-NPLTM 

IF(NOVT.CT.O)THEN 
DO  666  I-N1.N2 
666  I0UAD(I)“O 
ENDIF 

NZT-(NT0T»«2+NT0T)/2 

FHZ-FMC*1.0E6 

ETA-(376.7,0.) 

CAM-CMPLX(0. ,2 . «PI/WV) 

XK-2,»PI/UV 

q«.ooi«wv 

IF(NWR.EQ.O)A-Q 


SET  UP  SLEEVES,  NO  SLEEVES  (1,  0)  ON  THE  SECHENTS. 

DO  126  I-l.INM 
ISC(I)-0 
126  CONTINUE 

IF(NGO.EQ.0)THEN 

WRITE  OUT  WIRE/PLATE/OVERLAP  GEOMETRY  FOR  OKS  PLOTTING 

WRITE (9 , * ) NPLTS , NPLTH , NM , NP , NWR , NAT , WV , NOPL , NOVT 
DO  168  1«1,NP 
VRITE(fl,*)X(I) ,Y(I) ,2(1) 

168  CONTINUE 

DO  169  I-l.NM 
WRITE(9,»)IA(I) ,IB(I) 

169  CONTINUE 

DO  161  NPL-1, NPLTS 

161  WRITE(9,*)NCHRS(NPL) ,NPL11(NPL) ,NPL22(N?L) ,NDNPLT(NPL) ,IPN(NPL) 
DO  152  I-l.NPLTM 
DO  163  >1,4 

WR1TE(9,*)PA(I,J,1) ,PA(I , J ,2) ,PA(I , J ,3) ,PB(I , J , 1) ,PB(I ,  J,2) , 

2  PB(I, J,3) 

163  CONTINUE 
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162  continue: 

DO  166  NPL*1,NPLTS 
NCNR-NCNRS(NPL) 

DO  167  NC*1 ,NCNR 

WRm<9,*)PCN(l,NC,NPL)  , PCN (2  ,NC,NPL)  ,PCN(3,NC,NPL) 

157  CONTINUE 
166  CONTINUE 

DO  166  1*1 ,NOPL 

WRrr£(9,*>IOVT(I,l),IOVT(l,2)  ,I0VT(i  ,3)  ,I0VT(I,4)  ,ITK(I) 

166  CONTINUE 
END  IF 

C  END  OKS  PLOTTING  OUTPUT 

IF  <  NPRINT . NE . 0 . AND . NPRI NT , NE . 2 ) 00T06 1 0 
C 

C  WRITE  OUT  INPUT  DATA. 

C 

WRITE(6 , 62 1 ) FMC , WV , A , INTP , INTO , INT , IFIL 
621  FORMAT(/3X, ’OHIO  STATE  UNIVERSITY  ELECTROMAGNETIC  SURFACE  • 

2  ’PATCH  (ESP)  CODE:  VERSION  IV' 

2  //10X , ‘INPUT  DATA'//3X, 'FREQ. (MHZ)« ’ ,FlO . 3.3X, ’WAVE(M)= ' . 
2F10.3.3X, 'WIRE  RAD(M)*  ’.F9.6./, 

33X , ’ INTP* ’ i 14 , 3X , ’ INTO* 1 i 14 , 3X ( 1 INT  -  ',14 ,3X, ’IF!-*’ ,12) 
WRITE (6, 623)  CMM 

623  F0RMAT(3X,  'WIRE  CONDUCTIVITY  *  ,F6 . 2 .  *  MEGAMHOS/M’) 

WRITE (6 , 3006) IDWR , ITW2 , ICN , IPL , IAT , ITOT , ID2T , IDZTF , IDMI , 

2  IDZTI , ID2TFI , 1DMZI 

3006  FORMAT (/3X, 'SUMMARY  OF  ARRAY  DIMENSIONS '/3X, 'IDWR  *  ’,14/ 

2  3X,’ITW2  »  ' , I4/3X, 'ICN  -  ',I4/3X,’IPL  *  ’,14/ 

3  3X,’IAT  -  ' , I4/3X, 'ITOT  -  ' .I4/3X, 'IDZT  «  ',17/ 

4  3X, ’IDZTF  -  ' , I4/3X , ’ IDMI  -  ’, I4/3X , ‘ IDZTI  -  ’,17/ 

6  3X, ’ IDZTF I  -  ' , I4/3X, ’IDMZI  »  ’,12) 

WRITE (6 , 3006) IWRZM , IRDZM 

3006  F0RMAT(/3X, 'MATRIX  WRITE  AND  READ  Z  FROM  DISK  INDICATORS:*/ 
2  3X. 'IWRZM  -*, 12, 6X, 'IRDZM  -  ’,12) 

C 

C  WRITE  OUf  PLATE  GEOMETRY. 

C 

620  F0HMATC2X, 'X.Y.2  COOR .  (METERS)  OF  CORNER  ',12,'  -’.3F12.B) 


626  F0RMAT(//3X, 

'COOR. 

(METERS)  OF 

M3, 

,  ’  MODES  ON  PLATE  ’,13,/) 

626  F0RMAT(//3X, 

'THERE  ARE  ' 

.13.' 

MODES 

ON  PLATE  M3,/> 

630  FORMAT! 3X, 'I 

XAl 

YA] 

ZA1 

XA2 

YA2  ZA2  ’, 

2'  XA3  YA3 

ZA3 

XB1 

YB1 

ZB1 

XB2  ’. 

3 ’  YB2  ZB2 

XD3 

YB3 

ZB3 

V) 

660  FORMAT! IX, 13, IX, 18 (F4. 2)) 
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471 


474 


661 


662 

641 

472 

C 

C 

C 


k 

k 

k 


WRXT1(6,4?1)NPLTS 

FORMAT ( /3X , '  GEOMETRY  FOR  THE  '  .13,  •  PLATES’) 

IF(NPLTS.LE.G)G0T0610 

SO  472  NPL«1 .NPLTS 

IFdREC(NPL)  . EQ . 1 ) WRITE ( 6 , 4 7 3 ) NPL , N  CNRS ( NPL ) , S  EOM ( NPL ) 
TPN(NPL)  , IQS (NPL)  .ZSHT(NPL) 

F0RMAT(//26X, 'PLATE  NUMBER  ',13,'  (RECTANQULAR)  ’ 

/3X, ‘NUMBER  OF  CORNERS  ■  '.I3/3X, 'MAXIMUM  SEGMENT  SIZE 
2X, ' (WAVELENGTH)  ■  ' .F10.6/3X, 'POLARIZATION  INDICATOR' 
'GENERATIHO  SIDE  INDICATOR  ■  ',13/ 

3X, 'SHEET  IMPEDANCE  (OHMS/SQ.)  *  ',E11.4,**J  >,Eli.4) 
IFdREC(NPL)  , EQ . 0) WRITE ( 6 , 444) NPL , NCNHS ( NPL )  .SEOM(NPL) 
IPN(NPL) .IGS(NPL) .ZSHT(NPL) 

F0RMAT(///25X, 'PLATE  NUMBER  ',13,'  (POLYGONAL)' 

/3X. 'NUMBER  OF  CORNERS  «  ' .I3/3X, 'MAXIMUM  SEGMENT  SITE 
2X,' (WAVELENGTH)  -  ' .F10.S/3X, 'POLARIZATION  INDICATOR 
•GENERATING  SIDE  INDICATOR  «  ’,13/ 

3X, 'SHEET  IMPEDANCE  (OHMS/SQ.)  -  ',Ell.4,'+J  '.Eli. 4) 
DO  474  NCNR-l.NCNRS(NPL) 

WRITE (6 , 020) NCNR , PCN ( 1 , NCNR , NPL ) , PCX ( 2 , NCNR , NPL ) , PCN ( 3 
N2“NDNPLT(NPL) 

Nl-NPL-1 

NA-1 

IF(NPL.OT.l)  NA«NDNPLT(N1)+1 
IF(NPL.GT.l)  NN-N2-NA+1 
XF(NPL.EQ.l)  NN-N2 
IF(IWR.OE.  DTHEN 
WRITE  (6,626)NN,NPL 
ELSE 

WRITE  (e,626)NN,NPL 
ENDIF 
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IFdWR.LE.O)GO  TO  472 
WRm.(6,631) 

FDRMATdX,  'MONOPOLE  MODE  XI  Y1  Z1  X2', 

'  Y2  22  X3  Y3  Z3  X4‘, 

'  Y4  Z4  '/) 

DO  641  I-NA.N2 

WRITE(6,661)  I,((PAd,J,K),K«l,3),J*1.4) 

FORMAT  (4X, '  A  ' .3X.I3.2X, 12(F9. 6)) 

WRITE(6,Eb2)  I ,  ((PBd  ,  J  ,K)  ,K»1 ,3)  ,  J“1 ,4) 

FORMAT  (4X, ’  B  ' .3X.I3.2X , 12(F9. 6)/) 

CONTINUE 

CONTINUE 

WRITE  OUT  OVERLAP  MODES 

IF  (NOVT  .EQ.  0)  GOTO  610 
IMX-0 

DO  2039  I«1.N0PL 
IF(IOVT(I ,2)  .EQ.  0) 

WRITE(6 ,2031)ITX(I) ,IOVT(I , 1) , IOVT(I ,3) , IOVT(I ,4) 

IFdOVTd ,4)  .EQ.  0) 

WRITE(6 , 2031)  ITK(I)  , I0VT(I,3)  .IOVTd.1)  ,  I0VT(I,2) 
1F((I0VT(I,2)  .EQ.  0  .OR.  I0VTd,4)  .EQ.  0)  .AND, 

DOVL(I)  .LT.  0.1*WV)  WRITE ( 6 . 2032 ) DOVL (I ) 

IFdOVTd, 2)  .NE,  0  .AND.  IOVTd.4)  .NE.  0)THEN 
IF(IWR.LE.O)THEN 

WRITE(6,2033)ITK(I)  ,IOVT(1 ,1)  ,IOVTd  ,2)  ,I0VT(I,3)  ,IOVT(1 ,4) 
ELSEIFdWR.OE.  DTHEN 

WRITE(6 ,2030) ITK(I) , IOVT(1 , 1)  ,IOVT(1 ,2) ,  IQVTd ,3) , lOVTd ,4) 
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ENDIF 

ENDIF 

IF(IWR  .LE,  0  .OR.  ITK(X)  .EQ.  0)  GOTO  203? 

IMN»IMX+1 

IKX«IMX+ITK(I) 

DO  2040  ICT«IMN,IMX 
II-NPLTM-N0VT+1CT 

2041  IF(ICT.E(J,IMN)WR1TE(6,631) 

WRITERS, S61) XI,  ((PAUI.J.K)  ,K*1 ,3) , >1 ,4} 

VMTE(6, 662)11,  <(PB<n,J,K).K»l.  3)  ,>1.4) 

2040  CONTINUE 
2039  CONTINUE 

2030  F0RMAT(//3X,  ’CQOR,  (METERS)  OF  ',13.'  OVERUP  MODES', 
ft  '  BETWEEN '/3X, 'PLATE' ,13,',  SIDE' ,12, '  AND  PLATE ' ,13, 

A  ',  SIDE', 12/) 

2033  F0RMAT(//3X, 'THERE  ARE  ',13,'  OVERLAP  MODES’, 

ft  '  BETUEEN'/3X, 'PLATE', 13,',  SIDE' ,12,'  AND  PLATE', 13, 
ft  ',  SIDE',  12/) 

2031  FORMAT (//3X , 'COOR.  (METERS)  OF  ',13,'  OVERLAP  MODES', 

ft  '  BETWEEN  FACE  OF  PLATE ’.13.'  AND  PLATE  M3,’  SIDE  ',12/) 

2032  F0RMAT(3X ,'•*»»«  WARNING:  LENGTH  OF  OVERLAP  MODE  = 
ft  1PE10.2, '  LESS  THAN  0. 1 ‘WAVELENGTH ‘ ) 

610  CONTINUE 
C 

C  WRITE  OUT  LIST  OF  PLATE  SURFACE  PATCH  MODE  LOCATIONS 

C 

IF(IWR.GE.1)THEN 

VRITE(6,2046) 

2046  F0RMAT(//3X , ' LIST  OF  PLATES  ON  WHICH  SURFACE  PATCH  MODES  LIE’ 

2  /4X , ' N ’ , 3X , ' N-NWR ' , 3X , ' MON .  A',3X,'M0N.  B'/) 

D02046IL-1.NPLTM 

I«IL+NWR 

WRITE(6, 2047)1, IL.MPLA(IL) ,MPLB(IL) 

2047  FORMAT (IX, I 4 ,3X,I4,6X,I3,7X,13) 

2046  CONTINUE 

ENDIF 

IF(NPRINT.NE. 1 . AND . NPRINT . NE , 2)C0T0620 
IF (NWR . EQ . 0)G0TO726 
C 

C  WRITE  OUT  POINTS  AND  SEGMENTS. 

C 

WRITE (6 , 130)NP 

130  F0RMAT(///6X,I3, '  POINTS  ON  THE  WIRE  (METERS) ’ ,/SX , 'I ' ,9X , 

2 'X  (I)*,9X,*r  (I) ’ ,9X, *Z  (I)'./) 

D0140I-1.NP 

WRITE(6, 160)1, X«)  ,Y(I)  ,Z(I) 

160  F0RMAT(3X,I3,3(2X,E11 .4 ,2X)) 

140  CONTINUE 

160  FORMAT (///6X, 13,'  SEGMENTS  ON  THE  WIRE' ,/6X, 'J' ,4X, 'IA(J) ' , 

23X, ’IB(J) ' ,21, ’D(J)  (METERS)'/) 

WRITE (6 , 190) MAX , MIN , NWR 

190  F0RMAT(///6X, 'MODES  ON  THE  WIRE  STRUCTURE’ ,//6X, ‘MAXIMUM  NUMBER 
2  OF  MODES  AT  ONE  POINT  ■  ' ,12/61, 'MINIMUM  NUMBER  OF  NODES  AT  ONE 
30INT  -  ',12/BX, 'NUMBER  OF  WIRE  MODES  »  >,13/) 

IF ( IWR . LE . 0 ) G0T0729 
WRITE (6, 200) 

200  F0RMAT(/6X, 'I',2X, 'Il(I)’, 2X, >12(1) ’,2X, *13(1) ' ,2X, ' JA(I) > ,2X, 
a'JBd)',/) 

D0210I-1.NWR 

WRITS<6, 220)1,11(1),  12(1),  13(1)  ,JA(I),JB«) 


ESPO6620 

ESP06630 

ESP06640 

ESP06660 

ESP06660 

ESP06670 

ESPO66B0 

ESP06690 

ESP06600 

ESP06610 

ESP06620 

ESPO6630 

ESP06640 

ESP06660 

ESP06660 

ESP06670 

ESP06680 

ESP06690 

ESP06700 

ESP06710 

ESP06720 

ESP06730 

ESP06740 

ESP06760 

ESP06760 

ESP06770 

ESP067B0 

ESP06790 

ESP06800 

ESP06810 

ESP06820 

ESP06830 

ESP06840 

ESP06B60 

ESP06860 

ESP06870 

ESP06880 

ESP06890 

ESP06900 

ESP06910 

ESP06920 

ESP08930 

ESP06940 

ESP06960 

ESP06960 

ESP06970 

ESP06980 

ESP06990 

ESP07000 

ESP07010 

ESP07020 

ESP07030 

PESP07040 

ESP07060 

ESP07060 

ESP07070 

ESP07080 

ESP07090 

ESP07100 

ESP07110 
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220  FORMAT (2X ,6(13 ,41) ) 

210  CONTINUE 

720  WRITE(6, 160)NM 

230  F0RMAT(//6X, 'SEGMENT  LENGTHS (M) 1 ,/6X , ’ J ’ ,3X, 'IA(J) ' ,3X , ’ IB(J) ' > 
29X,‘D  (J)',/) 

DO  70  J-l.NM 
K-JA(J) 

L-IKJ) 

D(J)*SQRT((X(K)*X(L))»»2*(Y(K)-Y(L))»»2«(2(K)-Z(L))»*2) 

EGD“CEXP(GAM*D(J)) 

SSD ( J ) « (EGD- 1 ■ O/EGD) /2 , 0 
COD(J)«(EGD+1.0/EGD)/2.0 
WRITE(6 ,240)J , IA(J) , IB(J) ,D  (J) 

240  F0RMAT(3X,I3,4X,I3,6X,I3,3X,E13,6) 

70  CONTINUE 
C 

C  WRITE  OUT  ATTACHMENT  POINT  GEOMETRY. 

C 

726  1F(NAT .  LE .  OG0T062G 
WRITE (6, 660) NAT 

660  F0RMAT(/////3X,  "GEOMETRY  FOR  THE  M2,'  ATTACHMENT  POINTS1/) 
WRITE(6,570) 

670  FORMAT (3X,'l  SEGMENT  END  PLATE  B  (METERS)  V) 

D0B80NA-1.NAT 

NND-0 

IFCNSA(NA) .0T.NM)NND“i 
MMD-NND+1 

nas*nsa(na)-(mmd-i)»nm 

WRITE (6 , 690 ) NA ,  NAS , NND . NPLA (NA) , BOSK (NA ) 

690  FORMAT ( 2X , 12 , 3X , 14 , 6X , I 1 , 6X , 12 , 6X , F6 . 5 ) 

680  CONTINUE 
620  CONTINUE 

WRITE(6,722) 

722  F0RMAT(//6X, ‘LISTING  OF  LOADS  AND  GENERATORS'/) 

D0622J-1.NM 

IF(CAB3(VG( J)) . GT . 1 .OE-6) WRITE (6,723 )VG(J) , 3 
IF (CABS ( VG ( J ) ) . GT . 1 . OE-6 ) IANT- 1 
IF(CABS(2LD(J) )  .GT.l  .0E-6)WRITE(6,624)ZLDU) ,  J 
JJ*J*NM 

1F(CABS(VG( JJ)) .GT . 1 .0E-6)WR1TE(6,626)VG( JJ) , J 

IF ( CABS ( VG ( J J ) ) , GT . 1 . OE-6) IANT- 1 

IF  (CABS  (YU)  (JJ) ) .GT. 1 . OE-6) VRITE(6 ,626) ZLD(JJ) , J 

723  FORM.*! (3X , 2E13 . 4 ,  ’  VOLTS  BY  PT.  A  OF  SEGMENT  ',13) 

624  F0RMAT(3X,2E13.4, '  OHMS  BY  PT.  A  OF  SEGMENT  \I3) 

626  F0RMAT(3X,2EI3.4,  *  VOLTS  BY  PT.  B  OF  SEGMENT  M3) 

626  F0RMAT(3X,2E13.4, '  OHMS  BY  PT.  B  OF  SEGMENT  ’,13) 

622  CONTINUE 

D0627J-1.HAT 
IFCNAT .E0.O)C0T0627 

IF (CABS (VGA (3)) .GT. 1 .0E-6)WRITE(6,626)VGA(J) , J 

IF(CABS(VGA(J)) .GT.l. OE-6) IANT“1 

IF(CABS  (ZLDAO) ) .  GT .  1 .  OE-6)  WRITE  (6 , 629)  ZLDA  (J) ,  J 

628  F0RMAT(3X,2E13 .4,  ’  VOLTS  AT  ATTACHMENT  *,I2) 

629  FORMAT (3X,2E 13 . 4, '  OHMS  AT  ATTACHMENT  >,12) 

637  CONTINUE 

WRITE (6 , 833) NWR , NPLTM , NAT 

633  rORMAT(//6E,'NWR  -  NUMBER  OF  WIRE  MODES  •  ',14/BX, 

2  'NPLTM  ■  HUMBER  OF  PLATE  MODES  -  ',14/SX, 

3  'NAT  •  NUMBER  OF  ATTACHMENT  MODES  -  • ,14//) 

ICP-0 


ESP07150 

ESP07130 

ESP07140 

E8P071S0 

ESP07160 

ESP07170 

ESP07180 

ESP07190 

ESP07200 

ESP07210 

ESP07220 

ESP07230 

ESP07240 

ESP07260 

ESP07260 

ESP07270 

ESP07280 

ESP07290 

ESP07300 

ESP07310 

ESP07320 

ESPD7330 

ESP07340 

ESP07360 

ESP07360 

ESP07370 

ESP07380 

ESP07390 

ESP07400 

ESP07410 

E9P07420 

ESP07430 

E5P07440 

ESP07460 

ESP07460 

ESP07470 

ESP07480 

ESPO7490 

E3P07600 

ESP07610 

ESP07620 

ESP07530 

ESPO7640 

ESP07660 

ESP07660 

HSP07670 

ESP07660 

ESP07S90 

ESP0760C 

ESP07610 

ESP07620 

ESP07630 

ESP07640 

ESP07660 

ESP07660 

ESP07670 

ESP07680 

ESP07690 

ESP07700 

ESP07T10 
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IF (NFS1 , GT . 0 . AND . NFS2 . CT . 0) I CP* 1 

ESP07720 

IFUCP.EQ.OGO  TO  1226 

ESP07730 

1 

WRITER, 141) 

ESP07740 

141 

F0RMAT(//3X, ‘DEFINITION  OF  PORTS  FOR  MUTUAL  COUPLING  DATA'/) 

ESP077B0 

D0127NF-1.2 

ESP07760 

NFM-0 

ESP07770 

SN»1.0 

ESP07760 

NFS-NFSi 

ESP07790 

IF(NF.EQ.2)NFS-NFS2 

ESP07800 

D0128I-1.NWR 

ESP07B10 

IF(NFS.0T.NM)G0T0129 

ESP07820 

IF(I2(I) .NE, IA(NFS))GOT0128 

ESP07830 

NFM-I 

ESP07840 

.  * 

IF(JA(I).EQ.NFS)SN-1.0 

ESP078SO 

WRITE (6, 142) NF, NFS 

ESP07860 

142 

FORMAT (3X, 'PORT  ',12,'  IS  BY  PT.  A  OF  SEGMENT 

’,14) 

ESP07870 

G0T0128 

ESP07860 

0 

129 

IF ( 12 ( I ) . NE . IB (NFS-NM) ) C0T0128 

ESP07890 

NFM*I 

ESP07900 

IF(JBU)  . EQ . NFS-NM)SN«- 1 . 0 

ESP07910 

NFSNM-NF5-NM 

ESP07920 

WRITE (6 , 143) NF , NFSNM 

ESP07930 

143 

FORMAT (3X, 'PORT  ',12,'  IS  BY  PT.  B  OF  SEGMENT 

',14) 

E8P07940 

126 

CONTINUE 

ESP07960 

IF(NFM.GT.O. OR.NAT .EQ .0) GOTO 131 

ESP07960 

D0132I-1.NAT 

ESP07970 

IF (NFS . NE . NSA (I ) ) OCTC 132 

ESP07980 

NFM-NWR+NPLTM+I 

ESP07990 

WRITE (6, 144) NF, I 

ESP08000 

144 

FORMAT (3X, 'PORT  ',12,'  IS  BY  ATTACHMENT' ,13) 

ESP08010 

132 

CONTINUE 

ESP08020 

131 

CONTINUE 

ESP08030 

IF(NF.EQ.1)NFM1«NFM 

ESP08040 

IF(NF.EQ.2)NFK2-NFK 

ESP08060 

IF(NF.EQ.1)SN1«SN 

ESP08060 

IF  (NF .  EG .  2)  SN2-SN 

ESP08070 

127 

CONTINUE 

ESP08060 

1226 

CONTINUE 

ESP08090 

C 

ESP08100 

C 

WRITE  OUT  PARAMETERS  OF  FREQUENCY  SWEEP 

ESP08110 

C 

ESP08120 

IF(INDZI.NE.O)THEN 

ESP08130 

WRITE(6 , 3030) INDZI .FMC1 , FMC2 , DFZI , OFF 

ESP08140 

3030 

F0RMAT(/6X, 'PARAMETERS  OF  FREQUENCY  SWEEP  COMPUTATION’/ 

ESP081S0 

2 

3X, 'INTERPOLATION  METHOD!  INDZI  «  ',13/ 

ESP08160 

2 

3X, 'BEGINNING  FREQUENCY  (MHZ)  -  >,F9,3/ 

ESP08170 

3 

3X, 'ENDING  FREQUENCY  (MHZ)  -  >,F9.3/ 

ESP08180 

4 

3X, 'FREQUENCY  INTERVAL  FOR  COMPUTATION  OF  Z  (MHZ)  ■  '.F9.3/ 

ESP08190 

6 

3X, 'FREQUENCY  INTERVAL  FOR  MM  FIELD  COMPUTATIONS  (MHZ)  -  ' ,P9.3)ESP08200 

IF(IABS(IRS12)  .EQ .  1)WRITE(6 ,30B0) THRU, PHR1) 

ESP08210 

# 

3060 

FORMAT (3X. 'RADIATION  ANGLE:  THETA  »',F6.1,' 

PHI  «' ,F6. 1, '  DEO’ 

)ESP08220 

IF(IABS(IRS12).EQ.2)VRITE(8,30»0)THIHC.PHINC 

ESP08230 

3090 

F0RMATC3X, 'INCIDENT  ANGLE:  THETA  »',F6.1,'  PHI  »',F«.l,'  DEG’) 

ESP08240 

IF  (IABS  (IRS12) .  EQ .  2)  WRITE(6 ,  SIOOITWU) ,  PKRD 

ESP08260 

» 

3100 

FORMAT (3X, 'SCATTERING  ANGLE:  THETA  »',P6.1,’ 

PHI  -\F6.i, 

ESP08260 

*' 

DEG') 

ESP08270 

ENDIF 

ESP08280 

IF  (NGO .  Eq .  0)  G0T0600 

ESP06290 

NFZ-0 

ESP06300 

3010 

IF (INDZI .NE.O)THEN 

ESP08310 

59 

777 
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NFZ«NF2+1 

FHC-FMCl+(NFZ‘l)*DPZI 
WV*300 .O/FMC 
XK“2,0»PI/WV 
FH2-FMC*1.0E6 
OAM*CMPLX ( 0 , 0 , XK) 

Q»0,001»WV 
END  IF1 

COMPUTE  IMPEDANCE  MATRIX 

READ  IN  IMPEDANCE  MATRIX  IF  IRDZM  NOT  ZERO 
IF  ( IRDZM . LT . 1 ) G0T063 1 
IF(IFIL.EQ.0)READ(12)  (ZT(I)  ,I»1,NZT) 

IF(IFIL.EQ.1)THEN 
DO  777  J-l.NTOT 
DO  777  I-1.NT0T 
READ(12)ZTPd.J) 

END  IF 
CONTINUE 

IF (IRDZM. EQ • 3)CQT0276 

IF  A  FREQUENCY  SWEEP  IS  BEINC  MACE  SCALE  THE  1MAQINARY  PART  OF 
THE  PLATE  SHEET  ADMITTANCE  BY  THE  FREQUENCY 


636 


D0636NPL-1.NPLTS 
ZSHTF(NPL)-ZSHT(NPL) 

IFdNDZI.EQ. 0)0070636 

IF (ABS (AIMAO (ZSHT(NPL) )) ,LE . 1 , OE-34)C0T0636 
YSHT-l.O/ZSHTCNPL) 

FRATI0-FMC/FMC1 

YSHTF“CMPLX (REAL ( YSHT ) ,FRATIO*AIMAG(YSHT)) 

ZSHTF(NPL)-1.0/YSHTF 
CONTINUE 

SUBROUTINE  ZTOTZ:  COMPUTES  IMPEDANCE  MATRIX  ZT,  OR  ZTF,  AND  IS 
LOCATED  IN  ESP4SUBS  FORTRAN 

CALL  ZTOTZUA.IB.INM.ISC.Il.IQ.IS.JA.JB.MD.NVR.ND, 
2NM1NP1C(1D,SCD1D,X1Y,Z,ZLD,NPLTS,NAT,ZS,.TRDZM1ZLDA1 
3PA , PB , NSA , NPLA , PCN , IPL , IPLM , BDSK . ZT , ZTF , NM12N ,HM23N , ICN , 

4NDNPLT .NOVT , INT, INTP . INTO , CMM ,ERVSR , RMIN , DR , IAT . IPN , 

6IQUAD , NCNRS , IFIL , IREC , I CC , IERVSR , PDI ST , MPLA , MPLB , PSZ , ZSHTF , SEGMX ) 
276  CONTINUE 

IF  INDZI  •  3.  FORM  AND  STORE  EFFECTIVE  INVERSE  OF  Z 

IFdNDZI .  EQ .  3JTHEN 
CJ(i)»(l. 0.0.0) 

IFvIFIL.EQ.O)  CALL  SQROT(ZT,CJ ,0, 1 ,NTOT) 

CJ(1)-(1. 0,0.0) 

IF(IFIl.EQ.l)  CALL  CROVKZTF.CJ,  IDZTF,  1,0,1  ,NTOT) 

ENDIF 

IF(INDZI.NE.O)THEN 
IF (12468 . EQ . 0) IZC (NFZ) -NFZ 
IF (12468 , EQ . 2468) THEN 
003260X11-1,3 
IF(IZCdll)  .EQ.DTHEN 
IZHXT-III 
IZCdID-3 
ELSE 

xzcam-izcdiD'i 


ESP08320 

ESP06330 

ESP08340 

ESP08360 

ESP08360 

ESP08370 

ESP06360 

ESP08390 

ESP08400 

ESP08410 

ESP08420 

ESP08430 

ESP08440 

ESP0B460 

ESP08460 

ESP08470 

ESP08480 

ESP08490 

ESP08600 

ESP0B610 

ESP08S20 

ESP08B30 

ESP08B40 

ESP08660 

ESP08B60 

ESP08B70 

ESP08BB0 

ESP08B90 

ESP08600 

ESP08610 

ESP08620 

E5P08G30 

ESP0B640 

ESP08660 

ESP08660 

ESP08670 

ESP086B0 

ESP0B690 

ESPU8700 

ESP08710 

ESP08720 

ESP08730 

ESP06740 

ESP087S0 

ESP08760 

ESP0B77O 

ESP08780 

ESP08790 

ESP08800 

ESP08810 

ESP08820 

ESP08830 

ESP08840 

ESP088B0 

ESP0886C 

ESP08S70 

ESP08880 

ESP08BB0 

ESP08600 

ESP08910 
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ENDIF  EFP08920 

3260  CONTINUE  ESP08930 

ENDIF  ESP08940 

C  ESP06960 

C  IZCU)  ■  INTERPOLATING  POINT  NUMBER  OF  I  ENTRY  IN  ZT  OR  ZTF  ARRAESP06960 

C  ESP08970 

D03020J-I.NT0T  ESP08980 

IKdFIL.EQ.OUJ-J  ESP08990 

IP(IFIL.EQ.1))J«1  ESP09000 

D03O2OI-JJ.NTQT  ESP09010 

IF(IFIL.EQ.O)THEN  ESP09020 

IJN-I J(1 , J ,NTOT)  ESP09030 

IF (12468 ,  EQ . 0) ZTIN (NF2 , IJN ) -ZT  <1 JN)  ESP09040 

IF (12466. EQ. 246S)2TIN(IZNXT,IJN)-ZTdJN)  ESP09060 

END  IF  ESP09060 

IFdFIL.EQ.l. AND.  12468. EQ.0)ZTFIN(NFZ,1 , J)-ZTF(I , J)  ESP09070 

IFdFIL.EQ.l  .AND.  12468 1 EQ  ,2468)ZTFIN(IZNXT,1 , J)-ZTF(I , J)  ESP09080 

302C  CONTINUE  ESP09090 

ENDIF  ESP09100 

C  WRITE  OUT  IMPEDANCE  MATRIX.  ESPOSllO 

IF (IWRZM. LE. 0)G0T0632  ESP09120 

C  OPEN (UNIT- 1 ,NAME- ’ZMAT . DAT’ ,TY.  £= 'UNKNOWN ' , FORM-  UNFORMATTED')  ESP09130 

IF(IFIL.EQ.0)WRITE(12MZT(I),I«l,NZT)  ESP09140 

IFdFIL.EQ.UTHEN  ESP091S0 

DO  778  J-l.NTOT  ESP09160 

DO  778  I-l.NTOT  ESP09170 

77e  WRITE  (12)  ZTF.d  ( J)  E3P08180 

END  IF  ESP09190 

C  CLOSE (UNIT-1)  ESP09200 

632  CONTINUE  USP09210 

IF (IWRZT . GT . 0 . AND . IF IL , EQ , 0) WRITE(6 , 634) FMC  ESP09220 

634  FORMAT C//3X, 'LOWER  TRIANGULAR  PART  OF  SYMMETRIC  IMPEDANCE',  ESP09230 

2  '  i'.'TRIX  AT  ',F9.3,‘  (MHZ)  '//BE, 1 1  ‘,6X,  ’  J’  ,12X,  'Z(I,  J)  ’/)  ESP09240 

IF ( IWRZT . GT . 0 . AND . IFIL . EQ . 1 ) WRITE (6 , 636 ) FMC  ESP09260 

636  FORMA-d/m,  'IMPEDANCE  MATRIX  AT  \F9. 3,'  (MHZ)’//  ESP09280 

2  6X ,  ’ I  ’ , 6X ,  ,  12X ,  'Z(I ,  J) '/)  ESP09270 

D01234J-1 ,NTO 1  ESP092B0 

IF(IFIL.EQ.O) JJ-J  ESP09290 

1F(IFIL.EQ . 1 ) JJ-1  ESP09300 

D01234I-JJ ,NTOT  ESP09310 

IF  (IFIL. EQ.O) THEN  ESP09320 

IJN-IJ(I.J.NTOT)  ESP09330 

IF (IWRZT . GT . 0) WRITE( 6 , 1 233 ) 1 . J , ZT ( 1 JN )  ESP09340 

END  IF  ESP09300 

IFdFIL.EQ.l.  AND.  IWRZT.  GT.O.  AND.  J.EQ.J)WRITE(6, 1233)1,  J,ZTF(I,J)ESP09360 


1233  1 

F0RMAT(2I6,2E16.6) 

ESP09370 

1234  i 

CONTINUE 

E8P09360 

IF(I2466 . EO . 2468)GQT03260 

ESP09390 

IF (INDZI . NE . 0 . AND . NFZ . LT . 3) G0TD30 1 0 

ESP09400 

NFFS-0 

ESP09410 

NFFS-NFFS+1 

ESP09420 

3260 

CONTINUE 

ESP09430 

3060 

IF(INDZI.NE.O)THEN 

ESP09440 

IF(HFFS.EQ.1)THEN 

ESP0D4E0 

IF ( IRS 1 2 . Eq . 1 . AND , RF . LT . 0 . 0 ) TlffiN 

ESP00460 

WRITE(6,3110) 

BSP09470 

3110 

F0RMATC//3X, 'FREQUENCY  SWEEP  OF  ANTENNA  IMPEDANCE,  ' 

ESP09490 

2 

'EFFICIENCY,  AND  FAR- ZONE  GAIN'/ 

ESP0949O 

2 

2X,'F(MHZ) ' ,12X, 'ZIN(OHMS) ' ,8X,  "A  EFF' ,3X, 'GTHETA' , 2X, 

ESP09500 

3 

'  GPHI  ',21, 'GTHETA', 2X,'  GPHI  •) 

ESP09610 
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VRITE(6,3ti2) 

ESP09620 

3112 

F0RMAt(46X,  ‘•••HAS.  (D8)«<  *»PHASE  <DE0W/> 

ESP09B30 

ENDIF 

ESP09S40 

1FCIRS12  .EQ .  1 .  AND.  RF  .GE.  0.OTHEN 

ESP09550 

VRITE(6,3U3)RF 

ESP09660 

JUS 

F0RMAT<//3X, 'FREQUENCY  SWEEP  OF  ANTENNA  IMPEDANCE,  ' 

ESP0P570 

2  'EFFICIENCY,  AND  NEAR-ZONE  RAIN'/ 

ESP095B0 

3  31 ,  'PATTER.','  RADIUS,  RF  ■  \F11.6,'  METERS’/ 

ESP09S90 

1  2X,  ’F(MHZ) '  ,121,  ’ZIN(OKMS) '  ,BX,  '*/,  EFF’ ,3X, 'GTHETA' ,2X, 

ESP09600 

3  ’  OPHI  '  ,2X, '  3P.AD  ',2X,'0THETA',:X,'  OPHI  '.2X,'  ORAD  ’) 

ESP09610 

WRITEC6,3114) 

ESPCD620 

3114 

F0RMAT(46X,  ’••<  ••••MAG.  (DB) . .  ••••••PHASE  (DEG) . >/) 

ESP09630 

ENDIF 

ESP09640 

ENDIF 

ESP096B0 

c 

ESP09660 

C 

QUADRATIC  INTERPOLATE  ZTIN  OR  ZTFIN  IMPEDANCE  MATRIX  ARRAYS 

ESP09670 

c 

TO  FIND  THE  IMPEDANCE  MATRIX  AT  THIC  FMC 

ESP096B0 

c 

ESP09690 

112-1 

ESP09700 

c 

IF C INDZI  .EQ.:j)  112-2 

ESP09710 

DFF'i'=(i,F  'S-1UDFF 

ERP09720 

FMC-FMC1+DFFT 

ESP09730 

c 

TYPE*, NFrs, FMC 

ESP09740 

c 

RE-COMPUTE  F11EQ ,  DEPENDENT  QUANTITIES 

ESP09760 

W-300.0/FMC 

ESP09760 

XK-2.0*PI/VV 

ESP09770 

FHZ-FMC*1.0E6 

ESP09780 

OAM-CMPLX ( 0 , 0 , XK ) 

ESF09790 

Q-o.on:«Hv 

ESP09B00 

DO  3180  J-l.NM 

ESP09810 

IF CN  '  EQ. 0)03703180 

ESP09820 

K-UCJ) 

ESP09830 

L-1B(J) 

ESP09840 

D(J)«SQRT( (XOO-X(L))*«2  +  (YOO-Y(L))**2+(2(K)-Z(L))*»2) 

ESP098B0 

egd*cexp(oam*d(J)) 

ESP09860 

SOD  C  J)“  (EOD-l  .0/E0D)/2 , 0 

ESP0987C 

CGD(J)“(EGD+1 . O/EGD) /2 , 0 

ESP098B0 

3160 

CONTINUE 

ESP09890 

C 

ESP09900 

112*1. 60+DFFT/DFZI 

ESP09910 

1F(II2.LT. 2)112-2 

ESP09020 

1F(II2.CT.NFZ1-1)II2«NFZI-1 

ESP09930 

IF(II2.NE. II2LST)THEN 

ESP09940 

12468-2468 

ESPC09SO 

I12LST-II2 

ESP09960 

COT03010 

ESP09970 

ENDIF 

ESP09980 

111*112-1 

ESP09990 

II3-II2+1 

ESP10000 

D03270I1I*1 ,3 

ESF10010 

IFaZC(III).EQ.l)JJl-III 

ESP 1 0020 

IF(IZC(in).EQ,2)JJ2*m 

ESF 10030 

IF(IZC(1II) .EQ.S) jJ3“III 

ESP.0040 

3270 

CONTINUE 

ESP 10060 

FM1“FMC1+(II1-1> •DFZI 

ES-:no60 

m-FMCl+(II2-l)*DFZI 

ESP10070 

FM3-FMC1'(II3-15«DFZI 

ESP 10060 

CL1*ALOC(FMO 

ESP <0090 

CL2213— AL0G(FM2**2/(FM1*FK3) ) 

ESP10100 

CL2l“AL0G(FM2/FHl) 

ESP10110 
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WVl«300.0/PMl 

XK1»2.0*PI/WV1 

WV2-300.0/FM2 

XK2«2,0*PI/WV2 

W3*300.0/FM3 

XK3»2 , 0»P1/WV3 

DFZI2*DF21»»2 

EL1»0 . 6* (FMC-FM2) • (FMC-FM3) /DFZ12 

EL2— (FMC-FM1)  •  (FMC-FM3) /DFZI2 

EL3*0 ,6»(FMC-FM1)*(FMC-FM2)/DF2I2 

IFdVRZT  . GT.O. AND . IFIL.EQ .  0)WRITE(6,634)FMC 

IFdWftZT.  CT,  0,  AND ,  IFIL.EQ .  1) WRITE (8,636 )FMC 

WVMIN*300 .0/FMC2 

RIJ-0.0 

DQ3060J*1 iNTOT 

IFdNDZI.EQ.2)THEN 
C  FIND  CENTER  COOR.  OF  MODE  J 

IFU.LE.NURJTHEN 
I2J*I2(J) 

XCJ-XC2J) 

YCJ«Y(I2J) 

ZCJ*Z(I2J) 

ENDIF 

IF ( J . CT . NWR . AND , J . LE , NWR+NPLTM ) THEN 
JJJ*J-NWR 

XCJ«0.6*(PAUJJ,111)+PA<.)JJ,4,1))  . 
YCJ-0.S*1PA<JJJ,1,2)+PA(JJJ,4,2)) 

ZCJ-0 , 5» (PA  (JJ J  ,  1 ,3)  +PAQJ J  ,  4 ,3) ) 

ENDIF 

IF (J . CT . NWR+NPLTM) THEN 

JJJ-J-NPLTM-NWR 

JSB-NSA(JJJ) 

IFCJSG .LE.NM) I2J*IA(JSG) 

IF(JSG.GT.NM)I2J-JB(JSG-NM) 

XCJ*F(I2J) 

YCJ-YU2J) 

ZCJ-ZU2J) 

ENDIF 

ENDIF 

IFdFIL.EQ.0)JJ*J 
IFdFIL  EQ- 1)  JJ»1 
D03O60I-JJ.NT0T 
C 

C  IF  INDZI  *  2.  FACTOR  OUT  THE  EXP(-J*K*R)  DEPENDENCE  OF  THE 

C  ELEMENTS  IN  THE  IMPEDANCE  MATRIX  BEFORE  INTERPOLATING 

C 

CEXP1.(1. 0.0.0) 

CEXP2-U.  0,0.0) 

CEXP3-(1. 0.0.0) 

CEXPK-<1. 0,0.0) 

IFdNDZI  .Et).2)THEN 
C  FIND  CENTER  COOR.  OF  MODE  I 

IF(I.LE.NWR)THEN 
I2J-I2(I) 

XCI-X<I2J) 

YCI*Y(I2J) 

ZCI«Z(I2J) 

ENDIF 

IF (I . CT . NWR . AND . I . LE , NVR+NPLTM)THEN 
JJJ-I-NWR 


ESP10120 
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XCI*0.6*(M(JJJ,  1 , 1)+P*(JJJ  ,4,1)) 
YCX-0.6+<PA(JJJ1l12)+PAUJJ,412)) 
2CI*0.6+(PA(J.!.J,113>+PA(.!J.M13)) 

ENDIF 

IF(1  .QT.NWNPLTIDTHEN 

JJJ-I-NI’LTH-NVR 

JSG-NSAUJJ) 

If(JS0.U:.NM)I2>IA(J3G) 

IF(JSG ,GT .NM)I2J-IB(JSG-NM) 

xc:-xa2j> 

YCX-Y(I2J) 

ZCI-ZU2J) 

ENDIF 

RlJ-SQRT((XC.J-XCI)+*2+(YCJ-Ye:)»+2+(ZC.J-ZCI)**2) 

IF(RIJ/WHiN.GT .0 . 601)?HEN 
CEXP1»CEXP(-XJ+(XK1-XK2)+RIJ) 

CEXP2“CEXP(-XJ»(XK2-XK2)«RIJ) 

CEXP3“CEXP(-XJ+(XK3-XK2) *RIJ) 

CEXPK-CEXPI-XXXK-XKD+RIJ) 

ENDIF 

ENDIF 

IFUFIL,EQ.O)THEH 
XJN-XJd ,J|NT0T; 

FAC1-CEXP1 
FAC2-0EXP2  1 ' 

FAC3-CEXP3 

FACK-CEXPK 

ZTd.'N)"El.l+2TJN(  JJ1 , 1 JN)  /FAC1+EL2»ZTIN(J  J2  (IJN)/FAC2+ 

2  EL3*ZTIN  ( J  J3 , 1  Jf!)  /FAC3 
ZT(I JN)«ZT(IJN)*FACK 

IFCINDZI . EQ. 2. AH!), RIJ/WVhIN  LE .0 . 601  .AND . I , LE. KWR.1ND . J , LE.NWR) 
2  THEN 

XXXl"AIMAa(ZTIN( J J1 i IJN) ) 

XXX2-AIHAG(ZTXN(JJ2,IJN)) 

XXX3-AIMA0 (ZTIN U J3 , 1 JN) ) 

BBB-  (XXX3-2 . 0-XXX2+XXX1)  /CL2213 
CCC- (XXX2-XXX1 -BBB+CL21 ) /DFZI 
AAA-XXXl-BBB+CLl-CCC+FMl 

XXX- AAA+BBB* ALDO (FMC) +CCC+FMC 
2TdJN)-CMPLX(REAL\ZTdJN))  ,XXX) 

ENDIF 

IF(TVRZT.flT.0)VRITE(6, 1233) I, J ,ZT(IJN) 

END  IF 

IF«m.EC|.l)THEN 

FAC1-CEXP1 

FAC2-CEXP2 

FAC3-CEXP3 

FACK-CEXPK 

ZTF(I,  J)-EL1»ZTFIN(JJ1 ,1 ,  J)/FAC1+EL2»ZTFIN(  JJ2  |X  i  J)/FAC2+ 

2  EL3*ZTFIN(JJ3,I,J)/FAC3 
rrFCi.j)-ZTFa,j)*FACK 
XFdNDZl  .EQ.  2.  AND.  RIJ/WNXN.LF.  0.601)  THEN 
XF(IUnZT.CT.0)VRITE<6. 1233)1,  J.ZTFd.J) 

XXXI- AINAG (ZTFIN (J.)l , t  ,J) ) 

XXX2-AIMA0 (ZTFIN ( J J2 , I , J ) ) 

XXX3«AIMA0(2TFINCJJ3,I  >.!)) 

BBB- (XXX3-2 . 0+XXX2+XXX1 ) /CL2213 
CCC- (XXX2-XXX1-BBB+CL21 ) /DFZI 
AAA»XXX1-BBB*CL1-CCC»FM1 
XXX-AAA+BBB* ALOG (FHC) +CCC+FMC 
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2TF  ( I ,  J )  *CMPU  (REAL  { ZTF  (1 , 3  > ) ,  XXX)  ESP1 1320 

ENDIF  ESP II 330 

IF<mZT.aT,0)Witf?E(6,1233)i,1llZTFa1J)  BSP11340 

ENDIF  ESP113B0 

3080  CONTINUE  ESP11360 

ENDIF  ESP11370 
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COMPUTE  COUPLING  BETWEEN  TWO  WIRE  PORTS 

IF(ICP,£Q.l)  CALL  COUPLE (ZT , ZTF , NFM1 , NFM2 , SN 1 t SN2 f I 12 , V , 

2  NTOT.IFIL.ICC) 

IF(ISCAT.OE.i)  GOTO  B01 
COMPUTE  CONSTANT  VECTOR  A  SOLVE  SYSTEM 
IFdWR,GT.0)WRITE<6,676) 

678  P0RMATC//6X, 'ANTENNA  MODAL  CURRENTS ' /3X , 'MODE' ,3X, 

2  »REL  MAG',6X,'ABS  MAO' ,6X, 'PHASE' ,10X, '»**  COMPLEX  •*•'/) 

IF  (IANT  .Et) .  0)Q0T0600 

« 

FOR  ANTENNA  PROBLEMS ,  SET  UP  RMS  VECTOR  AND  SOLVE  FOR  CURRENTS 

CALL  ANTVdl , I2|I3, IA, IB ,IWR, JA, JB,NM,ZT, IFIL, 

A  ICC,  ZTF, CJP,  CG,  VG,  Yll,  ZU.VIN,  NGEN, NWR.NPL1M, NAT,  VGA, PIN, 

2  A.CMM.D.DISS.GAM.SOD.ZLD.ZS.ZLDA.INM.MD.ND.HSA.m) 

IF  THERE  IS  ONLY  ONE  GENERATOR  COMPUTE  INPUT  IMPEDANCE 

ZIN«(0. 0,0.0) 

YIN-(0. 0,0.0) 

IFCNGEN.EQ.DTHEN 

YIN’Y11/(CABS(VIN))*»2 

ZIN-1. O/YIN 

ENDIF 

112-2 

PRAD-PIN-DISS 

EFF-100.0-PRAD/PIN 

IF(INDZI.EQ.0)WRITE(6,B13)  Yll.2ll.EFF 

613  F0RMAT(/3X, 'INPUT  ADMITTANCE (MHOS)  -  '.F10.6,'  J  '.F10.6/ 

2  3X, 'INPUT  IMPEDANCE  (OHMS)  -  '.F10.3,'  J  '.F10.3/, 

3  3X,  'EFFICIENCY  (PERCENT)  -  \F7.3/) 

IF (INDZI.EQ. 0. AND .NGEN .EQ. 1) WRITE (6, 613)  YIN.ZIN.EFF 
IF (INDZI . EQ . 0 . AND . NGEN . NE . 1 ) WRITE (6 , 6 14 ) NGEN , EFF 

614  FORMAT (/3X.' INPUT  IMPEDANCE  AND  ADMITTANCE  ONLY  COMPUTED  IF', 

:i  '  WIRE  HAS  ONE  GENERATOR  */3X,  ‘NUMBER  OF  GENERATORS  -  >,14/ 

0  3X, 'EFFICIENCY  (PERCENT)  •  *,F7.3/) 

601  CONTINUE 

IFCIFF.NE. 1)  GOTO  668 

672  F0RMAT(//3X, 'ANTENNA  PROBLEM,  ISCAT  -  ',16) 

673  F0RMAT(//3X, 'BACKSCATTERING,  ISCAT  -  ',16) 

$74  F0RMAT(//3X, 'BISTATIC  SCATTERING,  ISCAT  «  ’,16) 

677  F0RMAT(//3X, 'FORWARD  SCATTERING,  ISCAT-  ’,16) 

XFdKD/.I .  EQ .  0)  THEN 
IKISCAT .  EQ .  0)  WRITE  (6 , 672)  ISCAT 
IFdSCAT .  EQ .  1 )  WRITE  (6 , 673)  ISCAT 
IF(ISCAT . EQ . 2) WRITE(6 , 674) ISCAT 
IF (ISCAT. EQ. 3) WRITE(8, 677) ISCAT 
ENDIF 

IFdSCAT. NE,0)G0  TO  601 

GET  INPUT  FOR  PATTERNS 
ANTENNA  PROBLEM  ISCAT-0 

860  DtM-PI/180. 

IF  (IFE  .NE.  1)  GOTO  801 

PERFORM  ELEVATION  PLANE  RADIATION  PLANE  PATTERN 
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PHI-PWE-DTH 

CC  NPTS*360/FNDFE*1.6 
CC  XNDFE-360/ (NPTS- 1 ) 

NPTS-ELRANa/FNDFE+1.5 

CC--SET  XNDFE»0,  INITIALLY  AND  CHECK  NPTS.CT.l  (AJF) 

XNDFE-0 , 

IF  ( NPTS , GT . 1 )  XNDFE-EIRANQ  /  (NPTS  - 1 ) 

IEA-1 

IF (NPLOTS .  CT .  0 ,  AND .  INCZI ,  EQ ,  0)  WUTE(8 ,338)  IEA ,  NPTS  ,PHFE 
338  FORMAT (IX, 2 (13, IX) , 2X,F7,2) 

IFCXNDZI ,NE,0)THEN 
PHI-PHRD-DTH 
NPTS-1 
END  IF 

DO  903  1*1 ,NPTS 

CC---ADDED  ELMIN  TO  THETA  IN  NEXT  LINE  TO  START  AT  ARB,  DBS,  PT, 
THETA-ELMIN  -DTH  ♦ < I - 1 ) -DTH- XNDFE 
IF (THETA . OT . 1 .0000 1 *PI )  THEN 
THETA-2. 0*PI-THETA 
IF( 19999 . NE . 9999)PHI-PHI»PI 
19999-9999 
ENDIF 

1FUNDZI .NE.O)THETA-THRD»DTH 

C  CALL  SORTS (IA, IB , II ,12 ,13,NWR,NM, A,CCD ,SQD,FH2,D , 

C  A  0,112, 1SCAT , ZTF , ZT , IFIL, ICC , ETT , EPP , 

C  A  X , Y ,Z ,NPLTS ,NAT ,PA ,PB ,NSA ,NPLA,PCN ,BDSK , IQUAD , 

C  A  NPLTM , IPL , IPLM , C JP , CJT , ETTS , EPPS , ETPS , EPTS , THETA , PHI , JA , JB , 

C  A  SC8P , SCST , SPPH , SPTM , STPM , STTH , IMAGE , ICN , NDNPLT) 

CALL  SORTBN (IA, IB, II, 12, 13, NWR , NM , A , COD , SOD , FHZ , D , 

•A  0, Il2,ISCAT, ZTF, ZT, IFIL, ICC, ETT, EPP. INTP.INTD, 

A  X , Y , Z , NPLTS ,NA7,PA,PB,NSA, NPLA . PCN , BDSK , IQUAD , 

A  NPLTM, IPL, IPLM, CJP.CJT.ETTS, EPPS, ETPS, EPTS, THETA, PHI, JA.JB, 

A  SCSP . SCST , SPPM , SPTM , STPM , STTM , IMAGE . ICN , NDNPLT , 

A  RF , EXN , EYN , EZN , EXT , EYT , EZT , EXP , EYP , EZP , ERRS , ETRS , EPRS , STRM , SPRM) 
PET ( I ) -PHS (ETTS ) 

PEP(I).PHS(EPPS) 

ETEd)-CAB3(ETTS)—2/<  <1.0, 0.0) -30.  O-PIN) 

EPE(I)-CABS(EPPS)--2/< (1.0, 0,0) -30. O-PIN) 

AET-AMP(ETE(D) 

AEP-AHP(EPEd)) 

AER-0.0 
PER(I)-0.0 
IF (RF . OT . 0 . 0) THEN 
PER (I) “PHS (ERRS) 

ERE(I)“CAB8(ERRS)«  2/(<l.  0,0, 0)*30.  O-PIN) 

AER-AMP ( ERE (I ) ) 

ENDIF 

1F(NPLOTS.GT.O. AND . INDZI . EQ . 0)URITE(8 ,333)DB(AET) , PET(I) , 

2  DB(AEP) , PEP ( I ) , DB( AER) , PER(I) 

II-I 

903  CONTINUE 

1F( INDZI .NE.O. AND.RF.LT. 0.0)TKEN 

HRITE(6,3120)FMC,ZIN,EFF ,DB(AET) ,DB(AEP) ,PET(II) ,PEP(II) 

WRITE (10, 312B)FMC,ZIN,EFF,DB(AET) ,DB(AEP) ,PET(II) ,PEP(II) 

3120  FORXAT(1X,F9.3,2X,E10.4, ’  J  • ,E10.4,2X,F8,1,2(2X,F6.2), 

2  2(2X,F6. 1)) 

3126  F0RMAT(1X,F9.3,2X,E10.4,2X,E10.4,2X,F8. 1 ,2(2X,F6.2) , 

2  2(2X,F6. 1) ) 

GOT03070 

ENDIF 
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IPdNDZI . NE. 0. AND. RF.CE. 0.0) THEN 
WRITE(6,31Si)FMC,ZIN,EFF,DB(AET) ,DB(A1P) ,DB(AER) , 

2  PET(n),PEPdi),PiRdi) 

WRIT£(lO, 3126)FMC,ZIN,EfF,DB(AET)  , DB(IEP)  ,DB(AER) , 

2  PET(XX) , PEPCII) ,PER(II) 

3121  FORHATax.FS.S^X.filOd,’  J  '  ,E10.4,2X,F6, 1 ,3(2X,F6.2) , 

2  3(2X,F6.1)) 

3126  FORMAT (IX  ,F9 . 3 , 2X ,E10. 4 , 2X.E10 .4 , !X ,F6 . 1 ,3(2X ,F6 . 2) , 

2  3(2X|F6.1)) 

ODTQ3070 
END  IF 

IF (RF . LE . 0 . 0) THEN 
WRITE (6, 904)  PHFE 

804  FORMAT!///. ’  FAR-ZONE  CAIN  ELEVATION  PLANE  PATTERN,  PHI  -’.F6.1. 
2  *  DEC. ’/) 

URITE(6,902) 

902  FORMAT!'  (DEC)  ••MAO  (DB)»»  *PHASE  (DEC)*'/ 

2  '  THETA  CTHETA  CPHI  CTHETA  CPHI'/> 

ELSE 

VRITE!6,804)  RF , PHFE 

804  FORMAT!/// , •  NEAR- ZONE  GAIN  ELEVATION  PLANE  PATTERN /3X , 

2  'R  •  ' ,F11 . 6, II , ’METERS  PHI  - > ,F6 . 1 . *  DEO.’/) 

URITEI6.80B) 

806  FORMAT!'  (DEO)  ••••••MAC  <DB)*»»«»*«*  •••••PHASE  (DEQ>»»*«*' 

2  ’  THETA  CTHETA  CPHI  OR  OTHETA  CPHI  OR’) 

ENDIP 

DO  906  I-l.NPTS 

CC<<<ADDED  ELM IN  IN  NEXT  LINE  TO  START  AT  ARB.  CBS.  PT. 
XII-ELMIN+(I-1)»XNDFE 
AET-AMP(ETEd)) 

AEP-AMP(EPEd)) 

AET-DB(AET) 

AEP-DB(AEP) 

IF!RF.CT.O.O)THEN 
AER“AMP(ERE(I) ) 

AER-DB(AER) 

URITE(6  ,B07)XIl ,  AET.  AEP .  AERiPETd )  ,PEP(I)  , PER(I) 

807  FORMATdX.FB.l.SdX.Fe^MISX.Fe.l)) 

ELSE 

URITE (6,907)  XII. AET .AEP, PET (I). PEP (I) 

907  FORHATdX.FB.l^X.Ffi.O.CX.Fe.S.SX.Fe.l.CX.Fe.l) 

EN01F 

ETE! I ) -CSQRT (1 . OE-20+ETE (I)) 

EPE  (I )  «-CSqRT  ( 1 .  0E-20+EPE  ( I ) ) 

IF(RF.CT.O.O)ERE(I)“CSt)RTd  .OE-20*ERE(I)) 

906  CONTINUE 
901  CONTINUE 

IF (IF A  .NE.  1)  GOTO  699 
C 

C  PERFORM  AZIMUTH  PLANE  RADIATION  PATTERN 

C 

THETA-tHFA*DTH 

C**NOTE:  IN  NEXT  TWO  LINES  CHANGED  360  TO  AZRANG 
NFTS-AZRANG/FNDFA+1 . s 
CC*«*SET  XNDFA-O.  AND  CHECK  NPTS.CT.l 
XNDFA-O. 

IF(NPTS . OT . 1 ) XNDFA“AZRANC/ (NPTS- 1 ) 

IEA-2 

IF(NPLOTS , OT . 0 , AND . INDZI . EQ , 0) VRITE(8 , 338) IEA , NPTS , TUFA 
DO  911  I®1, NPTS 
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CC***ADDED  AZMIN  IN  NEXT  LINE 

PHI»AZMIN*DTH  +(I-1)*DTH*XNDFA 

ESP13090 

ESP13100 

C  CALL  60RTB(IA,IB,11,I2,I3,NWR,NM,A,CGD,S0D,FHZ,D, 

ESP13110 

C  k  0,112, XSCAT , ZTF , ZT ,XF1L , ICC , ETT , EPF , 

E3P13120 

C  kX , Y , Z , NPLTS , NAT , PA , PB , NSA , NPLA , PCN , BOSK , 1QUAD , 

ESP13130 

C  kNPLTM , IPL , IPLM , CJP , CJT,ETTS , EPPS , ETPS , EPTS , THETA ,PHI , JA , JB , 

ESP13140 

C  k  SCSP , SCST , SPPM , SPTM , 3 TPM , STTM , IHAflE , 1 CN , NSNPLT) 

ESP13160 

CALL  S0RT8N (IA, IB, 11, 12, 13, NWR , NM , A , COS , SOD , FH2 , D , 

8SP13160 

k  0,H2,ISCAT|ZTF,ZT,IFIL,ICC,ETT ,EPP ,INTP ,INTD, 

ESP13170 

kX , Y , Z , NPLTS , NAT , PA , PB , NSA , NPLA , PCN , BOSK , IQUAD , 

ESP131B0 

kNPLTM , IPL , 1PLM , CJP , CJT , ETTS , EPPS , ETPS , EPTS , THETA , PHI , JA , JB . 

ESP13190 

A  SCSP , SCST , SPPM , SPTM , STPM , STTM , IMAQE , ICN , NDNPLT , 

ESP13200 

* 

k  RF , EXN ,  EYN , EZN , EXT , EY  T , EZT , EXP , ET P , EZP , ERRS , ETRS , EPRS , STRM , SPRM ) 

ESP13210 

PET ( I ) -PUS ( ETTS ) 

ESP13220 

PEP (I)»PHS( EPPS) 

ESP13230 

ETAZ( I)“CABS (ETTS) ••2/((1.0,0.0) *30.0*PIN) 

ESP13240 

• 

EPAZd)«CABSCEPPS)**2/(<1.0,0.0)»30.0»PIN) 

ESP132S0 

AET-AMP(ETAZd)) 

ESP13260 

AEP-AMP(EPAZd)) 

ESP13270 

AER-0 . 0 

ESP132B0 

PER«)-0.O 

ESP13290 

IF(RF.GT.O.O)THEN 

ESP13300 

PER(I) “PUS (ERRS) 

ESP13310 

ERAZ(I)*CABS(ERRS)»*2/ ( (1 . 0,0 .0) *30 . 0*PIN) 

ESP13320 

,  AER-AMP ( ERAZ ( I ) ) 

ESP13330 

END  IF 

ESP13340 

IF (NPLOTS . OT . 0 . AND . INDZI . EO . 0 ) WRITE ( 8 , 333 ) DB ( AET ) , PET (I ) , 

ESP13360 

2  DB(AEP) iPEP(I) ,DB(AER) ,PER(I) 

ESP13360 

911  CONTINUE 

ESP13370 

IF(RF.LE.OiO)THEN 

E5P13380 

WRITE(6,912)  THFA 

ESP13390 

912  F0RMAT(////2X,'FAR-Z0NE  GAIN  AZIMUTH  PLANE  PATTERN.  THETA  -> ,F6. 1ESP13400 

2, '  DEG.'/) 

ESP13410 

WRITE (6, 90S) 

ESP13420 

90S  FORMAT ( 1  (DEG)  »*MAG  (DB)**  »PHASE  (DEO)*'/ 

ESP13430 

2  1  PHI  GTHETA  CPHI  GTHETA  GPHI '/) 

ESP13440 

C  WRITEC6 ,914) 

ESP134S0 

C914  FORMAT ( '  PHI (DEG)  OTHETA(DB)  GPHI(DB) ’) 

ESP13460 

ELSE 

ESP13470 

WRITE(6,80B)  RF.THFA 

ESP13480 

BOB  FORMAT (/// , 1  NEAR-ZONE  GAIN  AZIMUTH  PLANE  PATTERN, >/3X, 

ESP13490 

2  'R  -  >  ,FU.6,1X, 'METERS  THETA  -'  ,F6 . 1 ,'  DEG. '/) 

ESP13600 

WRITE (6, 809) 

ESP13E10 

809  FORMAT ( '  (DEG)  ♦♦••••MAC  (DB)*»»**«»*  •♦•••PHASE  (DEC) ••••• '/  ESP13620 

2  1  PHI  GTHETA  GPHI  GR  GTHETA  GPHI  GR>) 

ESP13630 

ENDIF 

ESP13S40 

DO  913  I“1,NPTS 

XII“AZMIN+d*l)»XNDFA 

ESP135S0 

AET-AMP(ETAZ(I)) 

ESP13570 

AEP-AKP(EPAZ(I>) 

ESF13SB0 

AET-DB(AET) 

ESP 13690 

AEP-DB(AEP) 

ESP13S00 

IFCRF.GT.O.  OTHER 

ESP13610 

* 

AER“AMP(ERE(I)> 

ESP13B20 

AIR-DB(AER) 

ESP13630 

WRITE (6, 807)111, AET,AEP,AER,PETd)  ,PEP(I)  ,PER(I) 

E3P13640 

ELSE 

ESP136S0 

WRITE (6, 907)  XII,AET,AEP,PET(X) ,PEP(I) 

ESP 13660 

ENDIF 

ESP13670 

69 

o  o  o 


ETAZd)-CSQRT(l  , 0E-20+ETAZ ( I ) )  ESP13680 

EPAZd ) -CSQRT t 1 . 0E-20+EP AZ (I ) )  ESP 13890 

EftAZU )  -CSQRT  ( 1 . 0E-20+ERAZ(  1 ) )  ESP1 3700 

913  CONTINUE  ESP13710 

OOTD699  ESP13720 

601  CONTINUE  ESP13730 

ESP13740 

BACK  OR  BISTATIC  OR  FORWARD  SCATTER I NO  ISCAT  -  1  OR  2  OR  3.  ESP137B0 

ESP13760 

IFdNDZI  ,EQ.O)THEN  ESP13770 

IF(RF.U.0.0)WR1TE£6,721)  E8P13780 

721  PORMATOX.’FAR-ZONE  PATTERN')  ESP13790 

IF{RF.CT.O.O)WRITE(6,T31)RF  ESP13800 

731  FORMAT <3X, 'NEAR-ZONE  PATTERN:  R  ,F11. 6)  ESP13810 

IF ( IMAGE, EQ. 1) WRITE (6, 7 14)  ESP13820 

714  FORMAT (3X , 1  IMAGE  WAVE  INCLUDES ’ )  ESP13830 

IF( ISCAT ,EQ . 2) WRITE (6 ,7 13) THIN, PH IN  ESP13840 

713  FORMAT (3X, ’THETA  INC, (DEG.)  >  * ,F«, 1/3X, 'PHI  INC. (DEG.)  -  ’ ,F6. 1ESP13860 

2  //)  ESP13860 

ELSE  IF(NFFS.EQ.l)  THEN  ESP13870 

IF( IRS12 . EQ . -2) WRITE (6 ,714)  ESP13880 

WRITE(6 i 3200)  ESP13890 

3200  FORMAT (/ /3X , 'FREQUENCY  SWEEP  OF  TARGET  RCS  ’/  ESP13900 

2  12X, ’••*•••  MAO,  (DB/M**2)  *«•*••’,  E3P13910 

3  6X, ’••••••«  PHASE  (DEC)  . •»•’/  EgP13920 

4  2X , ’F(KHZ) ’  ,4X , ’STTM 1 ,4X, ’SPPM’ ,4X , 'STPM 1 ,4X, ’SPTM’ »  ESP13930 

5  6X i ’STTM 1 ,4X. ’SPPM’ ,4Xi ’STPM1 i4X, 'SPTM'/)  ESP13940 

END IF  ESP139B0 

ISP-0  ESP13960 

DTH-PI/180.0  ESP13970 

C  '  112*1  ESP139B0 

IEA-1  ESP 13990 

IF(ISE.EQ.O)C0T0916  ESP14000 

PHDG-PNSE  ESP14010 

PANG-PHDQ  ESP 14020 

NANG-O. 6+ (ELRANC/FNDSE)  ESP14030 

ANGRAN-ELRANG 

00T0917  ESP14040 

918  IEA-2  ESP140B0 

IF ( ISA . EQ . O)G0T0916  ESP14060 

THDO-THSA  ESP14070 

PANG-THDG  ESP14080 

NANG-0 . 8+ (AZRANG/FNDSA)  ESP14090 

ANCRAN-AZRANG 

917  CONTINUE  ESP14100 

C«*IN  LINE  BELOW  SUBSTITUTED  ANGRAN  (DEFINED  ABOVE)  FOR  360. 

CC***SET  DANG-O.  AND  CHECK  NANG.GT.O 
DANO-O. 

IF (NANG . GT . 0)DANG“ANGRAN/NANG  ESP14I10 

IA1-1  ESP 14 120 

ISP-ISP+1  ESP14130 

IF ( ISCAT • EQ  >  2  >  AND . ISP .  EQ , 1) 1 11  >0  ESP14140 

IA2-1+NANG  ESP14160 

NPTS-IA2  ESP14160 

IF(NPLOTS,GT.O. AND. INDZI.EQ.O) WRITE (8, 336) IEA,NPT2, ISCAT,  ESP14170 

2  pang. thin, phin  espuibo 

336  P0RMAT(1X,3(I3,1X),3(F7.2,2X))  ESPUIBO 

1F(INDZI . NE.O) IA2-1  ESP14200 


C»*NEW  LINES  TO  IMPLEMENT  CONSTANT  BISTATIC  ANGLE  OPTION  (TARGET  ROT.) 
C**IT  IS  ASSUMED  THAT  THE  REFERENCE  ANGLE  IS  THE  BISECTOR  OF  THE 
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* 


* 


C**BISTATIC  ANGLE.  NOTE  THAT  THE  SUM  OF  THE  INCIDENT  AND  OUTOOINS 
C**ANGLES  DIVIDED  BY  2  GIVES  THE  TARGET  NOTATION  ANGLE. 

DANCB-O, 

IF (NPTBIS .  OT .  1 )BANOB-BANGRG/  (NPTBIS-- 1 ) 

DO  1920  IBIS"1, NPTBIS 

IFUBISC  ,EQ .  1 )  THIN"  “BETA/2  .+DANOB*  (IBIS-1) 
1F(IB1SC.EQ.1)WRITE(S,92?8)IBIS,THIN 
9278  FORMAT! IX, 'IBIS* 1  (14, '  THIN-'  ,F12.3) 

DO  701  IANC"IA1|IA2 
CALL  GETCP2(ICPUB) 

CPU-ICPUB-ICPU 
WRITE (6 , 9SS2) 1AN0 , CPU 

9652  FORMAT (IX , ' XANG" ’ , I4,2X , 'CPU" ’ |F11 .2 |2X , 'SEC . ') 

CC* "ADDED  ELM IN .AZMIN  TO  THDG.PKDG  IN  NEXT  TWO  LINES 
IF ( IEA . EQ . 1 ) THDG-ELMIN* ( I ANC- 1 ) "DANC 
IFdEA.EQ.  2)  PHDC-AEMIN*  ( IANC- 1 )  "DANG 
IF( IBISC .EQ . 1 . AND . TEA . EQ . 1 )THDC-BETA/2 . +DANOB* ( IBIS- 1 > 
1FCIBISC.EQ.1. AND. IEA. EQ. 2)PHDG"B£TA/2 . +DANGB" (IBIS-1) 
THETA»THDfl»DTH 
PHI-PKDO*DTH 

IF(ISCAT.LE . 2) ISCT-ISCAT 
IF(IANG.EQ.0)ISCT*1 
IF( I ANO . EQ . 0)THETA-THIN"DTH 
IFCIANG . EQ. 0)PHI-PHIN*DTfi 
C  IFdWR.aT.O.AND.  1SCT.EQ.  1)URITE(6  ,  B76) 

III-l 

IFdSCAT.EQ  .3)  III"2 
DOTllIIS-l.III 
IF(ISCAT.LE.2)00T0716 
IFdlS.EQ.DTHEN 

IFCTHDG.GT- 180.0)THDOI"THDG-180.0 

IF(THDG . LE . 180 . 0) THDGI-THDG+180 . 0 

PHDGI-PHDG 

THETA-THDGI-DTH 

PHI-PHDGI-DTH 

ENDIF 

IFdlS  .EQ .  2)THEN 
THETA-THDG*DTH 
PHI-PHDG-DTH 
ENDIF 
ISCT-IIS 
716  CONTINUE 

IF(INDZI.NE.O)THEN 

IF(IANG.EQ.O)THEN 

THETA-THINC*DTH 

PHI«PHINC«DTH 

ENDIF 

IFdANG.EQ.DTHEN 

THETA-THRD*DTH 

PHI-PHRD"DTH 

ENDIF 

ENDIF 

C 

C  INSURE  THAT  (THETA, PHI)  IS  IN  THE  PROPER  RANOE 

C 

THETAP-THETA 

PHIP-PHI 

thdgo-thetap/dth 

PHDGO-PHIP/DTH 

IF(THETA . OT . 1 . 0000l«PI) THEN 


ESP14210 


ESP14220 

ESP14230 


ESP14240 

ESP14250 

ESP14260 

ESP14270 

ESP142B0 

ESP14290 

ESP14300 

ESP14310 

ESP14320 

ESP14330 

ESP14340 

ESP14350 

ESP14360 

ESP14370  . 

ESP 14380 

E3P14390 

ESP14400 

ESP144I0 
ESP14420 
ESP 14430 
ESP14440 
ESP14460 
ESP14460 
ESP14470 
ESP 14480 
ESP 14490 
ESP14600 
ESP14510 
ESP 14620 
ESP14E30 
ESP14640 
ESP14SS0 
ESP14660 
ESP14670 
ESP14680 
ESP14690 
E3P14600 
ESP 14610 
ESP14620 
ESP14630 
ESP 14640 
ESP14660 
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THETAP-2.0-PI-THETA 

PHIP-PH1+P1 

xf(phip , or . a . ooooi-pi ) phip-ph:p-2 .  o*pi 

THDOO-THETAP/DTH 

PHDOO-PHIP/DTH 

ENDIF 

C  TYPE1 1 1 1 .THETA/DTH , PHI /DTK , THETAP/DTH , PHIP/DTH 

1111  F0RHAT(1X,4F10.3) 

C  CALL  SORTS (1A, IB, II, 12, 13, HWR , HH , A , COD , SGD , FHZ , D , 

C  A  IWR.Itt.ISCT.ZTF.ZT.im.lCC.ETT.EPP. 

C  A  X ,  Y ,  Z , NPLTS ,  NAT , PA , P8 , NSA , NPLA . PCN , BDSK , IQUAO , 

C  A  NPLTM.IPL, IPLM.CJP ,CJT ,ETTS ,EPPS ,ETPS ,EPTS .THETAP ,PHIP , JA , JB , 

C  A  SCSP ,  SCST , SPPM ,  SPTM  , STPM ,  STTM .  IMAGE , ICN , NDNPLT ) 

CALL  SORTBN (1A, IB, II, 12, 13, HWR , HM , A , COO ,  SOD , FHZ , D , 

A  IWR, 112, ISCT.ZTF.ZT.IFIL, ICC, ETT  ,EPP,  INTP, INTO, 

A  X , Y , Z , NPLTS , NAT , PA , PB , MSA , NPLA , PCN , BDSK , IQUAD , 

A  NPLTM . IPL , IPLM , CJP , CJT , ETTS , EPPS , ETPS , EPTS , TKETAP , PHIP , JA , JB , 

A  SCSP . SCST , SPPM , SPTM , STPM . STTM . IMAGE , ICN , NDNPLT . 

A  RF , EXN , EYN , EZN , EXT , EYT , EZT , EXP , EYP , EZP , ERRS  p  ETRS , EPRS , STRM , SPRM) 
Til  CONTINUE 

IF(IANC.EQ.0)G0T0701 
710  FORMAT! IX , FlO .2.21, 2E20 .6) 

712  F0HMAT(//3X, 'SCATERR1NG  PARAMETERS*) 

17B  F0RMAT(2X,*  THETA- ’ .F10.2 ,  *  PHI- ' ,F10 . 2, '  SPPM-* ,E20,6, 

2*  sttm«*,e20.6,*  in  square-waves*) 

176  F0RHAT(2XI ’SPPM- • ,E20.6, *  DB  OVER  WAVELENOTH-SQUARED') 

173  F0RMAT(2X, 'STTM-' .E20.6, *  DB  OVER  WAVELENGTH-SQUARED*) 
PHTT-BTAN2( AIMAG (ETTS) .REAL (ETTS) ) -180.0/PI 
PHPP-BTAN2(AIMAG(EPPS) , REAL (EPPS) ) -180. O/PI 
PHTP-BTAN2 ( AIMAO (ETPS) .REAL (ETPS )) -180 . O/PI 
PHPT-BTAN2(AIMAG(EPTS) .REAL (EPTS) )• 180. O/PI 
PHTR-BTAN2 ( AIMAG (ETRS )  .REAL(ETRS) )  *180 .  O/PI 
PHPR-BTAN2 ( AIMAG (EPRS) ,REAL(EPRS))«180.0/PI 
IFCIANG.NEO.AND.NPLOTS ,GT.O. AND. INDZI .EQ.0)VR1TE(8 ,333) 

2  DB(STTM) .PHTT.DB(SPPM) .PHPP .DB(STPM) .PHTP.DB(SPTM) .PHPT, 

3  DB(STRM) ,PHTR,DB(SPRM) ,PHPR 

333  F0RMAT(6(1X,F6. 1 , IX , F6 . 1 ) ) 

ETAZ(IANG)-CMPLX(SQRT<STTM) . 1 .OE-20) 

EPAZ  ( I ANC)  -CHPL7.  (SORT  (SPPH) .  1 .  OE-20) 

STTM-DB(STTH) 

SPPM-DB(SPPM) 

STPM-DB(STPM) 

SPTM-DB(SPTM) 

STRM-DB(STRM) 

SPRM-DB(SPRM) 

IF (I ANG . EQ . 0 . AND . INDZI . HE . 0) GOT070 1 
IF( 1ANG . EQ . 1 . AND . IHDZ1 . NE . 0 ) THEN 

VRITE(6 . 3190) FHC . STTM , SPPM , STPM , SPTM .PHTT , PHPP , PHTP , PHPT 
3190  F0RMAT(1X,F8-2,4(2X,P6 .2) ,4(2X,F6. 1) ) 

IF (RF . LT . 0 . 0) WRITE ( 10 , 319$ ) FMC , STTM , SPPH , STPM , SPTM , 

2  PHTT. PHPP, PHTP. PHPT 
3196  F0HMAT(1X,F8.2,4(2X.F6.2)  ,4(2X,F6.  D) 

IF(RF.GE.0.0)WRITE(10,3196)FMC, STTM , SPPM , STPM . SPTM , STRM , SPRM , 

2  PHTT, PHPP. PHTP, PHPT, PHTR.PHPR 
3196  F0RHAT(lX,Fa.2,6(2X,F6.2)/9X,6(2X,F6.1)) 

GOT0701 

ENDIF 

IFCIANG.EQ. 1)THEN 
IFCRF.LE.O. 0>THEN 

write (6, 182) 


ESP14660 
ESP14670 
ESP14680 
ESP 14690 
ESP 147 00 
ESP14710 
ESP14720 
ESP14730 
ESP14740 
ESP 14760 
ESP14760 
ESP14770 
ESP14780 
ESP14790 
ESP 14800 
ESP14810 
ESP14B20 
ESP14830 
ESP14B40 
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ESP14860 
ESP14870 
ESP14880 
ESP  14890 
ESP 14900 
BSP14910 
ESP 14920 
ESP14930 
ESP14940 
ESP14960 
ESP14960 
ESP14970 
ESP14980 
ESP14990 
ESP1B000 
ESP1B010 
ESP1S020 
ESP15030 
ESP1B040 
ESP1B0B0 
ESP1S060 
ESP16070 
ESP1B080 
ESf’16090 
ESP16100 
ESP16110 
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162  FORMAT 1//2X ,  ’  *• (DEG)**  ••  CROSS  SECTION  (DB/H**2)  ♦* 

2  .  PHASE  (DEO)  *.•*•*>) 

write (6, 161) 

181  FORMAT(3X,'TH  PHI  STTM  SPPH  £TPM  SPTM  STTM', 

2  '  SPPM  STPM  SPTM 1 ) 

END  IF 

IP(RF.OT.O.O)THEN 

WRITE(6,179) 

ITS  F0RMAT(//3X, '•♦(DEO)**  •••••  NF.AR-ZOHE  CROSS  SECTION 

2  ' (DB/M**2)  *•*•*') 

WRITE(6, 1T8) 

178  FORMAT (3X, 'TH  PHI  STTM  SPPM  STPM  SPTM  STRM' , 

,  2  '  SPRM’) 

END  IF 
ENDIF 

IFdANC.OE.  1)THEN 

IF(RF.LE.O.O)WRITE(6,183)THDGO,PHDGO, STTM, SPPM, STPM, 

*  2  SPTM,PHTT,PHPP,PHTP,PHPT 

183  F0RMAT(1X,F6.1,1X,FB,1,4(2X,F6,2),4(2X,F6,1)5 
iF(RF.CT.O.O)URITE(6,184)THDGO,PHDGO, STTM, SPPM, STPM, 

2  SPTM. STRM, SPRM 

184  FORMAT ( IX, FS. 1,1X,FS,1,6(2X,F6,2)) 

ENDIF 

701  CONTINUE 
1920  CONTINUE 

C***END  OF  NEW  FIXED  BISTATIC  ANGLE  OPTION 
IF(INDZt.NE.0)G0T03O70 
916  CONTINUE 

IFdEA.EQ .  1)GOT0918 

699  CONTINUE 

CALL  GETCP2(KCPU) 

CC  CPU- (KCPU-JCPU)/ 100.0 

CPU-CKCPU- JCPU) 

WRITE (6 , SOS ) NRUN ,NWG , CPU 

608  F0RMAT(//3X, 'CPU  RUN  TIME  FOR  RUN '.IX, 13,'  GEOMETRY' , IX , 

2  13,'  «  ' ,F11.2, '  SECONDS'/) 

3070  CONTINUE 

IFdNDZI  .NE.O)THEN 
NFFS-NFFS+1 

IF (NFFS . It . NFF)G0T030S0 
ENDIF 

600  CONTINUE 

700  CONTINUE 

CALL  GETCP2(LCPU) 

CPU-(LCPU-ICPU) 

CC  CPU- (LCPU-ICPU) / 100 . 0 

WRITE (6. 607) CPU 

607  FORMAT (//3X, 'TOTAL  CPU  RUN  TIME  -  ',F11.2,*  SECONDS'/) 

9374  STOP 

»  END 

SUBROUTINE  GETCP(ICPU) 

C  CALL  TIMES (DATE, TIME, IVCPU.ITCPU) 

ICPU-ITCPU*0.0013 

*  RETURN 
END 

SUBROUTINE  CETCP2(ICPU) 

REAL  ETZME,TAJUUY(2) 

TIME-ETIME (TARRA Y ) 

ICPU-TARRAY(l) 

RETURN 


ESP 16260 
ESP16270 
ESP 16280 
ESP16290 
ESP 16300 
ESP16310 
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ESP 16330 
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ESP 16360 
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ESP1S370 
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ESP1S410 
FSP16420 
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ESP1B4S0 
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ESP16600 

ESP16610 

ESP16620 
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B8P16640 
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ESP1BE60 
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ESP1B680 
ESP16B90 
ESP16600 
ESP16610 
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ESP16660 
ESP16660 

ESP16670 
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ESP16760 
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END 


•  f  il«  (Unui  •»p4nM.r»trpl»t»  <untllt*d  plttt) 

ftRNCTftl 

NGO-1 .NPRINT-2 .NRUNS- 1 ,NV<33«1 , IUR-0 , IVR2T-0 , INT-4 , INTP-6 , INTD-1B , 
INVR-0 ,  IRGM-1 ,  IF1L-0  ,RF— 1 .0,  INDZI-0 , 

i 

kFSWEEP 

ft 

kpkrm 

IFE-0 , IPFE- 1 , FNCFE-3 , 0 , FHFE-90 , 0 , 

IFA-0 , XPFA-1 , FNDFA-3 , 0 .THFA-90 . 0 , 

I SE-2 , IPSE" 1 , FNDSE-3 , 0 , PHSE-0 .0,THIN"90.0,PHIN-0.0, EUUNG-0 . , 

ISA-0 , IPSA" 1 , FNDS A-3 . 0 , TMSA-90 . 0 , 

IBISC-1, BETA-120. «NPTBIS-121,BAN0RG— 360. , 

k 

AFVIRET 

FHC-1 300 . 0 , CMM-38 . 0 , A-0 . 00 1 , NPLTS- 1 , 

k 

APLATEC 

NCNRSU)-4,SEGMU)«0.2,IRECa)«0,IPN(l)*3,IGSU)-0,ZSHT(l)-(0. 0,0.0) 
XPCD  —  .4445,0. 4446,0. , 

YP(1)“0.  ,-.3429,0. ,.3429, 

ZP(1)*0. 0,0. 0,0. 0,0,0, 
ftPLATEO 
k 

ASAVEZ 

1WRZM-0.IRDZM-3, 

ft 

ftVIREAG 

NM-3 , NP-4 , NAT- 1 , NFPT- i , NFS 1 -0 .NFS2-0 , 

X(l) -0.0, 0,0, 0.0, -0.3, 

Y(l)-0. 0,0. 0,0. 0,0.0, 

Z(l)-0. 0.0.26,0. 6,0. 26, 

IA(1)-1 ,2,2, 

XB(1)-Si,3,4, 

ft 

ftCENLOD 

IFMM(1)-3,IABB(1)-0,VLGG(1)-(0. 0,0.0) ,ZLL(l)“(50. 0,0.0) , 

ft 

ftATTACH 

NASAT(1)-1 ,IABAT(1)“0,NPLA(1)«1 , 

VGA(1)-(1 . 0,0 .0) ,ZLDA(l)-(0, 0,0.0) ,BDSK(l)-0 . 4 , 

ft 

***"d»t»  f  11a  filtnuii  •»p4n«m.r»trtiltp  (tiltad  platd) 
ftRNCTlU. 

NGO- 1 , NPRINT-2 , NRUNS- 1 . NWGS- 1 ,  IWR-0 , IVRZT-0 , INT-4 , INTP-6 , INTD-1 8 , 
INUR-0 , IRGM-1 , IFIL-0 ,RF--1 ,0, INDZI-0, 

ft 

ftFSUEEP 

ft 

ftPATTRN 

IFE-0 , IPFE- 1 , FMDFE-3 . 0 , PHFE-90 . 0 , 

IFA-0 , IPFA- 1 , FNDFA-3 , 0 , THFA-90 . 0 , 

I8E-2 ,  IPSE- 1 ,  FNDSE-3 . 0 ,  PHSE-0 . 0 ,  THIN-90 , 0 ,  PHIN-0 . 0 ,  EUUNC-0 . , 
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